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A classical computation of gravitational bremsstrahlung in ultra-planckian collisions of massive 
point particles is presented in an arbitrary number d of toroidal or non-compact extra dimensions. 
Our method generalizes the post-linear formalism of General Relativity to the multidimensional 
case. The total emitted energy, as well as its angular and frequency distribution are discussed 
in detail. In terms of the gravitational radius rs of the collision energy, the impact parameter b 
and the Lorentz factor in the CM frame, the leading order radiation efficiency in the Lab frame 
is shown to be e ~ (r s /6) 3(d+1) 7cm for d = 0, 1 and e ~ (r s /6) 3(d+1) 7 c ^- 3 for d > 2, up to a 
known d-dependent coefficient and a ln7 cm factor for d = 2, while the characteristic frequency of 
the radiation is ui ~ 7/6. The contribution of the low frequency part of the radiation (soft gravitons) 
to the total radiated energy is shown to be negligible for all values of d. The domain of validity of 
the classical result is discussed. Finally, it is shown that within the region of validity of our approach 
the efficiency can obtain unnatural values greater than one, which is interpreted to mean that the 
peripheral ultra-planckian collisions should be strongly radiation damped. 
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1. INTRODUCTION 

The scenario of TeV-scale gravity with large extra dimensions [lj, jsj is currently being tested experimentally 
at the LHC [3] . In this class of higher-dimensional gravity theories the Planck mass can be of the order of a TeV, 
provided that the number of extra dimensions is greater than two. If so, the LHC will probe ultra-planckian 
physics where gravity not only is the dominant force [H , but in addition it is believed to be adequately described 
by the classical Einstein equations Q. This, presumably, allows one to make reliable theoretical predictions of 
gravitational effects without entering into the complications related to quantum gravity. One such prediction, 
namely the possibility of black hole formation at the LHC was extensively discussed in the literature during the 
past ten years [1, [6ri| . 

The main focus of this paper is the study of gravitational bremsstrahlung in ultra-planckian particle collisions. 
In ultra-planckian collisions with impact parameters smaller than the Schwarzschild radius of the center of mass 
energy of the colliding particles, the creation of a black hole is a dominant process and is accompanied by 
radiation losses at the level of 10-40% depending on the number of extra dimensions (see [13] generalizing 
earlier calculations by D'Eath and Payne [Til . Il2|). On the other hand, for impact parameters larger than 
the gravitational radius, the black hole formation process is exponentially suppressed and the gravitational 
bremsstrahlung becomes the dominant inelastic classical effect, which furthermore can be very strong due to 
the huge number of light Kaluza-Klein states in the graviton spectrum. It is worth noting, that gravitational 
radiation may not be maximal in head-on collisions. For instance, radiation from point particles falling into the 
black hole is minimal for radial infall and grows initially with the impact parameter, a fact which is well known 
both in four-dimensional general relativity and in models with extra dimensions. The way this is understood is 
the following: ultra-relativistic particles hitting with some critical impact parameter can be captured into quasi- 
circular orbits before they plunge into the hole. During these revolutions they radiate with almost constant rate 
and lose energy more efficiently than during a radial infall. 

A large number of papers from various groups is devoted to the study of gravitational bremsstrahlung in a 
variety of physical set-ups and using different approaches. A rather incomplete list includes the following: A 
classical calculation relevant to zero and small impact parameter using the colliding wave picture was suggested 
long ago by D'Eath [ll|,[l2j]. More recently this approach was extended to the extradimensional theories [9L . 
The approach of Smarr et al. in [l3j was based on the low frequency approximation for the amplitude to 
estimate the full energy loss. Matzner and Nutku proposed the use of the method of virtual gravitons [l4| in 
analogy with the Weizsacker- Williams approach in quantum electrodynamics. 

Another approach consists in considering radiation in the linearized theory from particles falling towards 
a black hole flH fl(| . Gravitational bremsstrahlung was also studied in the Born approximation of quantum 
gravity (in four dimensions ') 11711 . while the relationship between the Born approximation and the post-linear 
formalism was discussed in [18| . More recently, a novel approach, making use of perturbative quantum gravity 
to derive effective equations of motion and to describe gravitational radiation using those, was proposed under 
the name of "effective theory" approach in (l9j . Estimates of radiation losses in the eikonal approximation 
of (multidimensional) quantum gravity were given in while Amati, Ciafaloni and Veneziano in a series of 
papers [2(| have considered the problem of radiation in ultra-high energy collisions of massless particles in the 
context of quantum string theory. Estimates of the gravitational bremsstrahlung in high energy collisions in the 
ADD model were given in [2l[ . Finally, a lot of numerical work has also been devoted to this subject [2^ - [27| . 
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In four-dimensional theories, in particular, gravitational bremsstrahlung in small-angle ultrarelativistic scat- 
tering of point particles was calculated long ago, most notably by Kovacs and Thorne [28| . in the framework of 
the so-called "fast motion approximation" scheme, first proposed by Bertotti, Havas and Goldberg 29]. A simi- 
lar independent calculation, using momentum space perturbation theory up to second order in the gravitational 
constant, was performed in Ref. [l8[ leading to essentially the same results. The latter approach amounts to 
performing a perturbation expansion of the metric around the Minkowski background up to second order in the 
gravitational constant and is applicable for arbitrary mass-ratios of the colliding particles. In the case of one 
mass M much larger than the other m, the same result had been obtained earlier by Peters, also using a linear 
approximation but around the Schwarzschild background of the heavy mass [30| . 

The purpose of this paper is to study gravitational radiation losses in the case of ultra-relativistic collisions 
with large impact parameters, for which the scattering angle is small and calculations can be performed reliably 
within the post-linear approximation scheme of General Relativity. This paper continues a series of investiga- 



tions [31H3J dcvot :d to bremsstrahlung in flat space-time arising under non-gravitational scattering of charged 
particles 31 1 and scalar radiation in gravity- mediated particle collisions [32j | . Here we focus on the problem of 
gravitational radiation in collisions of massive point particles interacting gravitationally. The computation is 
purely classical and iterative and, as such, it can only be reliable within a certain domain of validity in the space 
of parameters. Somewhat unexpectedly, it turns out that, for parameter values within the region of validity of 
the classical approximation the radiation efficiency, i.e. the fraction of the initial energy which is radiated away, 
is substantially enhanced compared to the four-dimensional case, contrary to earlier qualitative estimates |34| . 
Moreover, for a number of extra dimensions greater than two, this quantity grows with a dimension-dependent 
power of the Lorentz-factor 7 of the collision. 

Specifically, for the leading order radiation efficiency in the laboratory frame we obtain (up to a dimension- 
dependent numerical coefficient, which however can be quite large, and an inessential In 7 factor for d — 2 extra 
dimensions) 1 

E /r< ? \ 3 ( rf + 1 ) 

7cm, «f=0,l (1.1) 



£ ~ V b J 



r„ \ 3{<J+1) 

1 %T 3 , d>2, (1.2) 



(?) 



where b is the impact parameter of the collision, r$ is the Schwarzschild radius for the center-of-mass collision 
energy 



1 



'TV 



d + 2 



(1.3) 



and £q = 771(7 + 1) — twy is the initial energy in the Lab frame. For collisions of two equal mass particles the 
Lorentz factors in the Lab and in the center-of-mass frames are related by 7 = 2^ m — 1. 

In the special case of d = our results for the total emitted energy as well as for its angular and frequency 
distribution are in perfect agreement with the results of Kovacs and Thorne [28[ . 

The present paper is organized in five sections of which this Introduction is the first. In Section 2 the 
model, the basic formulae and the iterative procedure, that will be employed, are described. Also, a complete 
orthonormal set of graviton polarization tensors in arbitrary dimensions, appropriate for the problem at hand are 
explicitly constructed. In Section 3 the expressions of the local and non-local contributions to the gravitational 
radiation source are obtained. The phenomenon of destructive interference of the leading local and non-local 
terms is shown. The total radiation amplitude is also obtained. Section 4 contains the computation of the 
leading ultra-relativistic order emitted energy in the classical theory in arbitrary dimensions. It also contains a 
detailed discussion of the frequency and angular distributions of the emitted radiation. The zero frequency limit 
of the radiation is analyzed in detail in all dimensions. In the same Section the domain of validity of the classical 
computation is discussed. The final result for the radiation efficiency is given and the unnatural possibility for 
it, to take values greater than one and even diverge in the massless limit for d 2 extra dimensions, is critically 
analyzed. A summary together with a few final comments are given in the final Conclusion section. 



1 Notice that these expressions differ from the generic formula e ~ (^s/fr) 3 ' d+1 ^7cm +1 given erroneously for all dimensions in [33ll . 
The difference is due to an error in [33| . related to the extent of the phenomenon of destructive interference discussed in the text. 
Nevertheless, the qualitative conclusions of l.'j.'J are still correct and will be discussed further in Section 5 of the present paper. 
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2. GENERAL SETTING 



We will study gravitational bremsstrahlung both in uncompactified Minkowski space-time M\,d-i = Mi^xM. 
and in the ADD model A/1,3 X T d with the d extra dimensions compactified on a torus. In both cases the collision 
will be confined on a brane M\$, and the D— dimensional cartesian coordinates are split as x M = (x^\y l ), 
x^ G Mi, 3, y l £ R d or T d , respectively. The brane subspace Mi i3 is selected by the initial conditions on the 
particle velocities and is fixed throughout the collision process, since the gravitational interaction cannot expel 
the particles from the brane. 



A. The model 



Consider two point masses m and m' moving along the world-lines x M — z m (t) and x M — z iM (t') (M = 
0, 1, . . . , D — 1) and interacting with the D— dimensional gravitational field. The corresponding action is 

S = -J2\ i (e 9M NZ M z N + —) dr-^ I R D J—gd D x, (2.1) 



where x|j = WitGd, e(r) is the ein-bein of the trajectory and the summation is over the two particles. The 
metric signature is chosen to be (+, — , — ) and our convention for the Riemann tensor is R B nrs = ^ns r ~ 
t nr,s + t ns T a r - r NR T AS> with T NR = ( l / 2 )g AB ' (9br,n + 9nb,r - 9nr,b)- Finally, the Ricci tensor and 
curvature scalar are Rmn = &a R A mbn and R = g MN Rmn, respectively. 

The action is invariant under the general coordinate reparametrization of spacetime, as well as under the 
independent reparametrizations of the particle trajectories. For the trajectory of particle m the transformation 
is r — > f = t(t), under which all fields arc scalars except for the einbein e(r), which transforms according to 
e(r) — > e(f) = e{r)df /dr. 

Varying S with respect to e(r) and z M (r) one obtains the equations 

e 2 g MN z M z N = m 2 (2.2) 

and 

d , , jy n e 
Tr {eZ 9MN ^ = 2 



ez N g MN ) = - g N p,MZ N z p , (2.3) 



respectively. Two analogous equations are obtained by varying with respect to e' and z' M , while variation of 
9m n leads to the Einstein equations 

G mn = 1 x 2 dT mn^ t mn = e f z M z N S D (^-z(r)) dT (2 4) 



Strictly speaking the notion of point-like particles is not compatible with full non-linear gravity. Nevertheless, 
it still makes sense in the context of the weak-field perturbation expansion around flat spacetime, which will be 
adopted here. In this approach one writes the metric as 

gMN = VMN + KDtlMN , (2.5) 

and expands all quantities in powers of HmNi using tjmn to raise and lower the indices. The flat background will 
be either D-dimensional Minkowski space M\ d—i> or the product Mi ; 3 x T d (d = D — 4) of four-dimensional 
Minkowski and a d— torus. 

It is convenient to define the quantity 

i>MN = h M N - 7}hri MN , h = ?] MN h MN , (2.6) 

and fix the general coordinate reparametrization symmetry by choosing to work in the flat-space harmonic 
gauge 

d N ^ MN = . (2.7) 

The reparametrization freedom of the particle trajectories will be dealt with later. The expansion of several 
relevant quantities in powers of Hmn is given in Appendix [X] In particular, the Einstein tensor in the harmonic 
gauge becomes 

yc >f 2 

Gmn = — 7£~ a i>MN ^-Smn + N MNl (2.8) 
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where □ = 9m<9 a/ is the flat d'Alembert operator, Smn is the 0(h 2 ) part of Gmn given by 

SmnW =hjtf (flNQ,P ~ hjip t Q) + h PQ (h,MP,NQ + h>NP,MQ ~ hpQ^MN — h.MN,PQ)~ 

~\ hPQ M h PQ,N - \h MN Uh + ^vmn (^2h PQ nh PQ - h PQ ^ L h PL < Q + ~h P Q, L h PQ > L j , (2.9) 

while Nmn stands for all cubic and higher in Hmn terms in the Einstein tensor. 

A few remarks are in order. Usually in the ADD scenario one assumes that matter is localized on the brane, 

1. e. restricted to the subspace M 13 . This assumption is non-contradictory within the linearized gravity, since the 
matter energy momentum in this case is purely flat-space tensor, whose conservation does not involve gravity. 
In the full non-linear gravity the Z?-dimensional Bianchi identity 

T A % = -i== Pm V—d\,N ~ \ 9 np,mT np = (2.10) 

implies D-dimensional geodesic equations for the particles (12.31) (and similarly for brane fields), which are 
generally inconsistent with the assumption of matter localization on the brane. For this reason we have to 
consider the world-lines in (I2.1[) from the beginning as D-dimensional curves. In general, without any extra 
forces ensuring confinement, the particle motion can be localized on the brane only by suitable initial conditions. 
However, it will be shown that the gravitational interaction treated perturbatively will not expel the colliding 
particles from the brane at least to leading order in our iterative scheme. This will be sufficient for the consistency 
of our computation of bremsstrahlung. 

Another technical difference from the standard ADD approach, also related to the fact that our treatment 
appeals to the .D-dimensional picture rather than to the four-dimensional one, will be the treatment of the 
graviton as a Z?-dimensional massless particle, instead of as a set of four-dimensional massive fields of spins 

2, 1, by suitable rearrangement of the components of the Z?-dimensional graviton Q. In particular, it will be 
more convenient here to choose directly the D{D — 3)/2 polarization tensors of the graviton propagating in the 
bulk, rather than to split them into polarizations of massive particles in four dimensions. 



B. The iteration scheme 

Our approach is identical to the one used in simpler models in [3ll | and (32j | . The iterative solution of the above 
field equations and gauge fixing conditions amounts to formally expanding and computing all fields step-by-step 
in the form 

$= °$+ 1 $ + 2 $ + ... , (2.11) 

where <& denotes any of the fields z M (r), z' M (t'), e(r), e'(r') and Hmn{x) and with the left superscript labeling 
the step of the iteration. 

The zeroth order solution is trivial. It describes the two particles moving with constant velocities (u = 
(u^, 0, 0) and u' M = (m' m , 0, 0), respectively in the Mi. 3 subspace, whose indices are labeled by lower case 
greek letters fi, v — 0, 1, 2, 3. Correspondingly, their trajectories °z M — ( °z a '(t), 0..., 0), and similarly for °z' M , 
are the straight lines 

V(r) = «"(0) + u^t , V V) = ^(0) + u'V, (2.12) 
in the absence of any gravitational field 

°h MN = 0. (2.13) 

The Lagrange multipliers are chosen equal to the corresponding particle masses 

°e = m, °e'=m', (2.14) 

so that the trajectories are parametrized by the corresponding proper times and the four-velocities satisfy the 
normalization conditions t]mnu m u n = u 2 = u' 2 = 1. 

We choose to work in the Lorentz frame in which the target particle m' is initially at rest and with the 
projectile m moving along the z-axis. Also, we introduce the space-like vector b = (6 M , 0, 0) lying on the 
brane, with: 

= ^(0)-z ,A1 (0). (2.15) 
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With no loss of generality one can further assume (bu) = (bv!) — and choose the x-axis along b. Thus, 

u» =7(l,0,0,u), u" 1 = (1,0,0,0) , 6 M = (0,&,0,0) , (2.16) 

where 7 = — v 2 and b is the impact parameter. 

The first order correction is obtained next. The zeroth order straight particle trajectories are sources of the 
first order gravitational field 1 Hmn{ 1 h' MN ) of each particle, which in turn, causes the first order deviation of 

the trajectory 1 z' M ( 1 z M ) of the other one. In the process, the first correction 1 e{ 1 e') of the einbein fields is 
also obtained. Explicitly, from the zeroth order trajectories one obtains the zeroth order energy-momentum 
tensor 

°r- = £ m /„ W(l -M^, (2.17, 

which in this order has only {M, N} = {/i, v\ components, and from the first order Einstein equations, given 
in the harmonic gauge by 

□ V MJV = -* D °T MN , (2.18) 

the first order correction 1 if) MN to the metric is obtained. Although only the 1 ip^ v components of ip MN are 
non- vanishing, the corresponding 1 h MN will also have non- vanishing bulk components from the trace part in 
Eq. (I2~18l) . 

Using 1 h MN and the zeroth order solution in equations (|2.2p and (12.31) one obtains for 1 e and x z M the 
equations 2 

l e = -y {hb 1 h M NU M u N + 2 VM nu m 1 z N ) (2.19) 



and 



( 1 eu M + m 1 z M ) = -x D mH P QM u p u Q , H PQM = 1 h PM ,Q - ^ 1 h PQM , (2.20) 



which upon elimination of 1 e give for z 

rrMN d i _ n MN tt ,,p,,Q to oi\ 

11 — Zn — —XdII tlpQN u , (.^-^tj 
dr 

where 

IL MN = V MN - u M u N (2.22) 

is the projector onto the subspace orthogonal to u M . Equations (I2.18|) . (|2.19p and (|2.2ip . together with two 
similar equations for 1 e' and 1 z' M , form a complete set of equations in this order. The gauge fixing condition 
(|2.7|) for 1 iIj mn is a consequence of the conservation of °T MN . 

It should be kept in mind that in the equation of motion (|2.21[) of particle m the gravitational field on the 
right hand side is the one due to ml . Singularities due to the action on m of its own gravitational field can be 
removed in lowest order by classical renormalization of the affine parameter of its trajectory |4(ij |: this will be 
enough for our purposes. In what follows we will omit self-action terms and consider only mutual gravitational 
interaction. 

In the next order of iteration one obtains the leading contribution to gravitational radiation. It is due to the 
accelerated particles, as well as to the cubic gravitational self-interaction. Indeed, the Einstein field equation 
reduces for the second order correction 2 ipMN of the gravitational field to 

□ 2 4>mn = -xd T M N , (2.23) 

with the source term having three contributions: 

tmn = 1 Tmn + 1 T' MN + Smn ( h) ■ (2.24) 



2 Our gauge condition is gMN% M z N = 1- To this order it reduces to x e = 0. 
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The first two terms of tmn (|2.24l) represent the particles' contribution to radiation. As obtained to this order 
from (|2.4p they are 3 



1 T M n{x) = m 



2 1 z {M u N ) + k d ^2u p 1 h P{M u N) y UmUn ^J ~ u mUn 1 z p d f 



5 D {x-°z(T))dT. (2.25) 



Ditto for m! with u replaced by it'. To repeat, in (|2.25|) 1 h is due to to'. Finally, the part Smn of T mn represents 
the contribution of the gravitational field itself to gravitational radiation. Given that Smn is quadratic in Hmn, 
consistency of the iteration to this order requires to substitute Hmn — > 1 h-MN inside S. Furthermore, just as the 
self-interaction terms were ignored in the particles' equations of motion, one has to keep in Smn( 1 h) only the 
products of the first order corrections h due to different particles. The various contributions to the gravitational 
radiation to this order is shown schematically in Figure [TJ 

It is straightforward to verify on the basis of the equations of motion satisfied by the first order fields that 

d N r MN = , (2.26) 

which guarantees the validity of the gauge fixing condition (|2.7[) to this order. 

So Smn considered as a quadratic form in Hmn constitutes the non-local (in terms of the flat space picture) 
source of gravitational radiation. This non-locality is due to the non-linearity of the theory and, as will be 
shown, leads to important differences in the radiation spectrum from linear theories like electromagnetism. 
More detailed discussions of this point in four dimensions can be found in [18], as well as in its generalization 
to arbitrary dimension but in a simpler non-linear model in (32|. 



1. The asymptotic behavior of the fields 

The same tensor Smn considered as a quadratic form in 2 Iimn plays the role of an effective energy-momentum 
tensor of gravitational waves [37| . This interpretation is based on the conservation equation (|2.26j) in combi- 
nation with the asymptotic properties of the solutions of the flat-space d'Alembert equation and is valid in all 
space-time dimensions. 

Recall, that in four dimensions 1 h liv , which according to (|2.18|) is due of a static or a uniformly moving mass, 
falls-off as 1/r, and the corresponding "field strength" (derivatives 1 h fJiL , y \) behave as 1/r 2 . On the other hand, 
the fall-off of the radiative component of the retarded potential 2 /i M „ is 1/r and the same for its derivatives 
2 h fl v t \. Correspondingly, in D-dimensional Minkowski space-time the behaviors of 1 1imn and its derivatives 
1 Hmn,l are l/r D ~ 3 and l/r D ~ 2 , respectively, where r is the radial coordinate in the D — 1-dimensional space. 
The corresponding fall-off for the radiation field can be found from the recurrence relations for Green's functions 
in neighboring dimensions. The derivation looks different in even and odd dimensions, but it leads to the same 
result, namely (see Appendix |D|) 

2 tlMN ~ 2 hMN,P ~ (£)_2)/2 ' (2.27) 

In the case of the ADD model the behavior of the field of a uniformly moving mass depends on the ratio of r 
to the compactification radius R. For r -c R it has the D-dimcnsional behavior l/r D ~ 3 , while for r 3> R it 



3 Symmetrization over two indices is defined according to A( MN ^ = (Amn + ^jvm)/^. 
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decreases as 1/r (with r being the radial coordinate in three-space). We assume that the observational wave- 
zone corresponds to this second condition. Then, the fall-off of the radiation field in the wave-zone will be the 
same as in four dimensions. 

Since the effective stress-tensor S , mat( 2 ^) is a homogeneous second order function of 2 h and its derivatives, 
its fall-off in D-dimensional Minkowski will be 1/A D ~ 2 ^ and in the ADD case 1/r 2 . This is precisely what is 
needed to get non-zero flux of radiation in the wave zone through the D — 2-dimcnsional sphere at spatial infinity 
of Mx d—u or through the corresponding two-dimensional sphere in the ADD scenario. Note that the use of 
pseudotensors to define the momentum density of gravitational radiation leads to the same result. Also, note 
that as a consequence of the above behavior the higher order non-linear terms contained in Nmn of Eq. (l2.8p 
do not contribute to the radiation energy-momentum. 



C. The energy-momentum of the gravitational radiation - Formulae 

Following Weinberg 37 1, we introduce the gravitational energy-momentum tMN = Gmn — G M ' N . Combined 
with the matter part Tmn, the total energy-momentum of the system is conserved, i.e. it satisfies 

d N {T MN +t MN ) = 0, (2.28) 

which is a consequence of the identity 8nG^ mn = 0, and which to second order in the above iteration scheme 
coincides with (|2.26[) . 



1. The uncompactified case Mi^-i 

Consider the fiat-space world-tube W bounded by two space- like hypersurfaces S±oo defined at t — > ±oo, 
and the infinitely far away time-like Cylindrical surface C . Thus, the boundary dW of W is the union of £-oo, 
S +oc and C. 

Using the energy-momentum conservation equation (|2 .28[) one can write successively for the emitted radiation 
momentum APm 



AP m = - f T MN da N + J T MN da N = - J T MN da N = - J t 

S + =c 2-^ dW W 

N+ jD„ _ I . i TV _ / . iN 



d N T MN d u x 



d N t M Nd u x= / t M Nda N = t MN da N . (2.29) 



w 



To justify the above steps, notice that the first equality is (with da > in both E_oo and S+oo surfaces) the 
statement that APm is minus the change of momentum of the matter system; for the second one uses the fact 
that the matter system is localized inside the volume and the contribution of C on the surface integral vanishes. 
Gauss' theorem was used next, together with (12.28)) to end up with the surface integral of tMN- Finally, for the 
last equality one used the fact that the surface integrals of tjfjv on E_oo and S+oo are zero. 

In particular, to the order of our computation in the previous subsection, using Gauss' theorem, equations 
and the gauge fixing condition, one may write instead of (|2.29[) 



AP M = J S MN ( 2 h)da N = J d N S MN { 2 h)d D x= l - J 2 h PQM a 2 ^ p Qd D x + f 2 h p Qu 2 ^ MP ^ Q d D x. 
c w w w 

(2.30) 

Furthermore, it is straightforward to show that the last integral on the right vanishes. Indeed, 

J 2 h PQ U 2 ^ M p.Q d D x = J 2 h PQ □ Vmp da Q - J 2 h PQ Q □ Vmp d D , 
w c w 

1 -J 2 h- p U 2 ij MP d D x 
w 

-i J 2 h □ 2 Vmf da p + X - J 2 h □ 2 ij M p' P d D x = 0, (2.31) 



c w 



where the first integrals of the first and third lines vanish as a consequence of the asymptotic behavior of the 
fields given in dD~4]) . (|D~5|) . and ([2T27]) . 
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Thus, the emitted momentum is 



AP M = \j 2 h PQM a 2 iJ; PQ dPx . (2.32) 
w 

Using the Fourier transformed quantities 

2 h MN (x) = tAzt J 2 h MN (k)e- ik * xQ d D k, ijMN ^ x ~j _ J T MN^ e -ik Q x^ d D k (2 33) 

one rewrites (|2.32l) in the form 

AP M - -^0y5 J k M G TCt {k)T SN {k)rl R {k)K SNLR d D k, (2.34) 

where rl {N (k) — tmn { — k) by the reality of tmn(x), and the tensor A SNLR is given by 

frSNLR = 1 ^SL^NR + ^R^L] _ _l_^Lfl (2 35) 

The retarded Green's function is G ro t = — V{\jk 2 ) + i7re(fc°)<5(fc 2 ). Its real part leads to an integrand, which 
is odd under parity k M — > —k and does not contribute to the integral. Thus, one writes equivalently 

AP M = 2(2 ^ D -i / 0(k )k M T SN (k)T* LR (k)A SNLR 6(k 2 ) d D k . (2.36) 

Taking into account the transversality of tmn (k M TMN{k) — 0) and the on-shell condition k 2 — of the emitted 
wave, one can replace the Minkowski metric in A SNLR by 

a m n — nM ttLN „mn i k M k N — 2(kg)k( M g N "> / Q7 \ 

= 9 L fc' 11 -V + j^p , (2-37) 

with any time-like unit vector g, where g n = 1 — g (g> g and fe/II = 1 + k' ® k' /(kg) 2 are projectors onto 
subspaces transverse to 3 and fe' = g IIfe = fc — (kg)g, respectively. Since (fc'g) = 0, the projectors ff II and 
fc'II commute. Their product A. M is then a symmetric projector onto the subspace M^ g , perpendicular to k 
and g. By construction, the projector A is idempotent (A 2 = A), thus on Mfc i9 it acts as the unit operator. In 
what follows, we will conveniently choose gu = u' M and calculate the flux in the Lorentz frame (the Lab) with 
< / = (l,0,...,0). 

Thus, A can equivalently be replaced in ()2.36|) by 

A SNLR = 1 [ A SL A NR + ^SR^L^ _ _}_ A SN A LR , (2.38) 

and then upon integration over |k|, one ends up with 



AP M = 4(2 ^ D _! J u D - z dio J dnk M T SN (k)T* LR (k)A SNLR , (2.39) 

S D-2 



where w = k° = |k|. It is easy to show that A = t}mn^- M = D — 2. Furthermore, A satisfies the following 
relations: 

(a) It is traceless on both pairs of indices 

VSN A SNLR = 0, VLR A SNLR = 0, (2.40) 

and (b) it is a projection operator, since it is idempotent (A 2 = A) 

7 ]P sriQTA MNST A PQLR = A MNLR , (2.41) 

and projects any second rank symmetric tensor, element of a space we denote by K, onto a subspace K\. Acting 
on a t mn (or t lr ) 4 , it returns a new symmetric tensor T ' MN = A mnp ®t P q with two indices. A acts as a 



4 The arguments below apply independently to the real and the imaginary parts of the complex- valued tmn that enters in 1 12. 391 1. 
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linear operator on the space K of symmetric tensors of rank two 5 . The properties of r' determine the space 
K A . 

In particular, the dimensionality of K\ is determined as follows: t' mn is a symmetric D x D matrix and 
thus it has D(D + l)/2 elements. In addition it satisfies the independent conditions i~' MN u' N = 0, r' MN k N = 
and rj MN t' mn = 0, which impose 6 D + (D — 1) + 1 = 2D conditions, leaving a total of D(D — 3)/2 indepen- 
dent elements in r'. Thus, the dimensionality of K\ is D{D — 3)/2, the same as the number of independent 
polarizations of the graviton in D— dimensions. 

Given that A acts on ATa as the unit operator one can write it (decomposition of unity) in terms of an 
orthonormal basis tensors {e^ N }, with V = 1,2, . . . , D(D — 3)/2 of Ka 

A MNLR = J24 IN 4 R - (2-42) 
v 

Using this expression for A, equation (|2.39[) takes the form 

AP M = 4(2 ^ g -i E / u D ~ 3d " J k M \r LR (k) , (2.43) 

v o s D - 2 

i.e. a sum of independent contributions one for each polarization. 



2. The compactified case Mi,3 x T" 



Following the same steps as in M\ £>_i , one is led to a similar expression for the emitted momentum in the 
ADD scenario. Consider the flat-space world-tube W bounded by two space-like hypersurfaces £±oo defined at 
t ±oo, and the infinitely distant time-like surface SV-j-ooj defined for — oo < t < oo, and whose intersection 
with the space-like hypersurfaces t — const will be denoted by B t . According to our convention all quantities 
constructed from the metric perturbation Iimn are flat-space tensors, so the covariant derivatives acting on 
them are simply partial derivatives. Starting with the energy-momentum conservation (|2.28p . and following the 
same steps as above, one concludes, first of all, that the change of the momentum of the system carried away 
by gravitational waves during the collision can be computed by integrating the flux of momentum through B t 
and also over time 



AP M = I dt J S MN (h)d D - 2 B N . (2.44) 



Note that the hypersurface B t is the product S 2 x T d and the corresponding measure is d D 2 Bm = r 2 dQ,2 d d y, 
with r the radial coordinate in R 3 . Thus, using Gauss' theorem one may write: 

oo 

J 8 N S MN d 4 xd d y= J S M0 d 3 xd d y~ J S M0 d 3 x d d y + hm / dt j ' d d y I S Mr r 2 dn 2 . (2.45) 

W Soo £_oo -oo T d 

Given that for t = ±oo the metric is flat 7 , the contribution of the integrals over S±oo is zero. Taking, in 
addition, into account the asymptotic behavior of 2 h and its derivatives (see Appendix [D]) , one concludes that 
([104]) takes the form: 

AP M = J d N S MN ( 2 h) d D x . (2.4(, ) 

w 



Finally, repeating the same steps and arguments as between equations (|2.30j) and (|2.32j) and taking into account 
the asymptotic behavior relevant to the case at hand, one ends up with the same expression as in the M% n—i 



5 This becomes even more evident if one denotes the pair of indices (MN) as a and the pair (LR) as b. Then the action of A on r 
takes the form r' a = A ab r b . 

6 The condition T~' MN k N = gives D — 1 constraints. This is due to the fact that k 2 = 0, which allows an arbitrariness in r' of 
the form akM^N- 

7 In this calculation we ignore the self-action of fields upon the particles, but even if we did not, the change of (infinite) self-field 
momenta would be zero. 
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Figure 2: The angles in lab frame used in the text. 



case for the emitted momentum 



AP M = i J 2 h PQM □ 2 i> PQ d D x . (2.47) 



w 

In terms of the Fourier transformed quantities defined by 



^ ' n 

(X ' V) = J2^V ^ / . (2.48) 



(V is the volume of the d— dimensional torus) the emitted momentum in the transverse to the brane directions 
takes the form 

Apl = I ~J dk ° dk ' E m^m^^wf ' (2 - 49) 

where k is a wave 3- vector on the brane and k% = n 2 /R 2 . AP l vanishes since the integrand is odd under 
parity n — > — n. The emitted momentum tangential to the brane can be obtained along the same lines as above 
with the only difference of using the discrete Fourier transformation on the torus. Finally, the emitted energy 
(E = APa) is 

OO 

E =Sv^J k2d ^J dn ™^ rl R (k) A— I , (2.50) 



where TMN(k) is the four-dimensional Fourier-transform of the total source evaluated at k° = y/k 2 + k\. 

The wave vector of the emitted gravitational wave is a genuine D-dimensional vector with quantized transverse 
components. We denote 

k M = jfe*,), k l T = —, n'eZ. (2.51) 

R 

The notation and conventions for the angles used in this paper is shown in Figure [2 For the brane four- vector 
W 1 we use the following parametrization ensuring zero square of k M : 

fc M = (^J vd 2 + k 2 ^, ms'mO cos ip, msmOsmip, tucos^^ . (2.52) 

Finally, using (|2.42l) the formula for the emitted energy becomes 

OO 

£= Jv£ s /■>**>/* i^w^r ■ (2 - 53) 
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D. Graviton polarizations 



The computation of the emitted energy can be simplified considerably if one chooses appropriately the or- 
thonormal basis of the D(D — 3)/2 graviton polarization tensors {s MN } that appear in (|2.43p and (|2.53[) . The 
construction of such a convenient basis is presented here. 

• The first step towards their construction is to define D — 2 space-like unit vectors orthogonal to k and u' 
and among themselves. The first D — 4 vectors e M , a = 3, . . . , D — 2 we choose to be orthogonal to the collision 
three-dimensional subspace spanned by the vectors u M , u' M , b M which lie on the brane and together with k 
fully characterize the kinematics of the problem under study. Thus, they are taken to satisfy the orthogonality 
properties 

e^ep M = -8 a p, (e a k) = {e a u') = (e a u) = (e a b) = . (2-54) 

As mentioned above, in the lab frame u' M — (1, 0, ... , 0). Also, with no loss of generality the vectors e a can be 
chosen to have vanishing components along the brane. 

Consider next that the brane Cartesian coordinate frame {t, x, y, z} together with the above D — 4 space-like 
unit vectors e^f form a D-dimensional right-handed basis {t, x, y, z, ea, . . . , e£>- 2 } and choose the £)-dimensional 
Levi-Civita symbol so that £ 0xyz3...(D-2) _ ^ The two remaining space-like vectors e\ and e 2 are defined as 
follows: e 2 is chosen orthogonal to both u M , u' M , namely 

ef = N- 1 e MMlMa '" MD - 1 u Ml u' Ma k Ma e !lMt . . . e (D _ 2) Md _, , (2.55) 

where N is the normalization factor which can be obtained by squaring the above expression. Specifically, 
making use of the relation k 2 = one obtains 

iV 2 = - [{u'k)u - {uk)u'f . (2.56) 

According to this definition 

(eau) = (e 2 u') = (e 2 fe) = (e 2 e a ) = , a = 3, 4, . . . , D - 2 . (2.57) 

Finally, the remaining vector e\ is chosen orthogonal to u' and to all previously constructed polarization vectors. 
It is given by 



ef = N- 



(ku) u' M - (ku') u M + (uu 1 ) - >- M 

[ku') J 



(2.58) 



with the same normalization factor. Note that e 2 is confined on the brane, while ef 1 has both brane and bulk 
components due to k M . 

• The next step is to use the above basis of orthonormal polarization vectors to build the N-p — D(D — 3)/2 
traceless mutually orthogonal polarization tensors e^ N , i.e. satisfying: 

ev N VMN=0 and e% N e v . mn = 5 V v> ■ (2.59) 
First, we build from e\ and e 2 two tensors, the ones familiar from four-dimensional gravity: 

p M p N _ p M p N p M p N , p M p N 

or, equivalently, the corresponding chiral ones 

-MN _ (I) (II) _ MN M _ e l ± le 2 f0 R1 \ 

e± -j= — e ± e ± , e ± — -j= . l^-O-iJ 

The remaining set of tensors is taken to have the same general form as the type-I and type-II above, namely to 
have (up to coefficients) the form of a sum of same index binomials ef 1 ef for the first kind, and different index 

binomials e[ M e^ for tensors of the 2nd type. Thus, the third tensor involving e\, e 2 , contains the diagonal 
sum of the D — 4 vectors e M 



E 



MN 



2(D - 4) 



(™5 ~ V D-2 



^Z? pM p N , M N 

M n e i e i t- e 2 e 2 
e„ e, 



D - 4 ^ 01 01 

Q = 3 



(2.62) 
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For D > 4 one can build the larger set of (D — 2)(D — 3)/2 tensors similar to erm, starting with the full set of 
unit vectors e\, e2, e Q , namely 



efef + efef 



J = 1,2,. 



Finally, the remaining D — 5 polarization tensors e Q are obtained from the linearly independent tensors 



MN 



1 



V(D~5)(D-A) 



M„N ST P M N 

&r e, 



{D-5)e%e»- ^ 



3,. 



(2.63) 



(2.64) 



built from e a , after ortho-normalization to satisfy the conditions (|2 . 59[) 8 . 
polarizations. 



For D C 5 there are no such 



• Finally, one can verify that \ MNLR can be written in terms of these D(D — 3)/2 polarization tensors 
e-p, V= I, II, III, (ij), ol, in the form 



k MNLR = Y, £ v N 4 R - 

V 



(2.65) 



As advertised, the advantage of using the above special choice of polarization tensors in (|2.43j) and (|2.53|) is now 
evident. The indices of the energy-momentum tensor tmn are carried by the four independent vectors it, it', &, k, 
of which the first three belong to the brane, while k M is multidimensional. These four entirely characterize the 
bremsstrahlung process. A consequence of the fact that the vectors e a are by construction orthogonal to those, 
is that only the three polarizations, namely with V = ± and III contribute to (|2.43p and (|2.53l) . Moreover, the 
e a terms in em can also be dropped and the latter simplifies to the "effective" polarisation: 



e (III) - 



D-A 



[efef 



2(D 



,M„iVl 



D-2 + - 



(2.66) 



Incidentally, notice that this effective vanishes, as it should, in the special case of four dimensions. 



Having chosen the vectors ef and e^ to be orthogonal to k and u , only the scalar products of ef and 
with u M and b will enter in the computation of the radiated energy momentum. According to (|2.57p . 



(e2ii) vanishes. The remaining non- vanishing products (ei&), (eitt) and (e2&) are given in the Appendix IA 2 



3. THE GRAVITATIONAL RADIATION SOURCE 



To linear order the metric perturbation is just the sum Hmn + h' RIN of the separate contributions of the 
two particles. They are for each particle the solution of (I2.18[) with the corresponding zeroth order energy 
momentum tensor (|2.17|) for its source. In Fourier space they are given by 9 



l-MN 



h MN(<i) 



2ir?CDm 

q 2 -Qt 



l " b 8{qu) 



u M u N 



h-MN — h' MN (q) 



h' MN {q) 



2-KXpvn! , , , 

„2 _ n 2 °\q u ) ( U MU N 

q Qrp 



d + 2 
1 



Vmn 



d + 2 



Vmn 



(3.1) 
(3.2) 



where q\ — l l / R, the transverse components of the momentum, and q\- = X)(9t) 2 - 

The deflection of the trajectory of particle m due to its gravitational interaction with m! is obtained by 



inserting the above h! MN in (|2.21[) and solving for 



1,M 



One obtains 



{2iryV 



v 2^ / a q i n 2 _ , 



(q 2 -q^){qu) 



JU 



iM 



,M 



d + 2 



7* _ qM 



2(qu) 



(3-3) 



8 One could a priori take a = 3, 4, . . . , D — 2 for the label of the polarizations f a . However, they are not independent, since as one 

D-2 

can easily verify y ' f£f = 0. 

a = 3 

9 It is worth noting that in four dimensions and in the center of mass frame the limit of /im jv f° r m ~ * an£ i 7 — > oo with £ = rwy 
fixed, is the Aichelburg-Sexl metric |3SH . 
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where 7* = 7 2 — l/(d + 2) . Similarly for 1 z' . As expected the trajectory deviation in the transverse directions 

<V~V / d*q / T 1 ^ qU [} e-^ b (e~^ s - l) (3.4) 

f^ d J r-Qr M v 1 



vanishes, since the summand is odd under parity (q l T —> — <7t0> an d as expected the particles do not leave the 
scattering plane. 



A. The local part 1 T M n + 1 T' MN of the source 



The contributions of the particle trajectories to the gravitational radiation source is given in (12.25[) . They 
can equivalently be written in the form 



l 2W(k) = e ikz ^m / dre tiku)T 



and 



1 T' MN {k)=e ikz '^m' / dre l{ku ' )T 



i(k 1 z) u M u N + 2U( M 1 z N ) - yt D —u M u N + 2x D u s ti s{M u N) 



i(k 1 z') u' M u' N + 2u[ M 1 z' N) - x D ~u' M u' N + 2x D u' s h s {M u' N) 



(3.5) 



(3.6) 



respectively. With the specific choice of polarization tensors made above, the projection of 1 T' MN (k) on any of 
them vanishes. Then, substituting the h! and 1 z given in (13.21) and (13.31) . one obtains: 10 



1 ?MJv(fc) 



,. ik b 



{2n) 2 V 



E 



' ku 2(ku) 2 



1 \ l*i 
(k ■ I ) ] u M u N + TT^-M(M-fjv) 



(ku) 



(3.7) 



in terms of the integrals I and r M , defined by 
rt _ [ S(pu') 8(ku - pu) 



P 2 -Pt 



,4 T i f S(pu')S(ku-pu)e-^ 4 
d P , Im= 9 5 Pm d p, 



P 2 -Pt 



with p l T = l l /R. These were computed in [42| and expressed in terms of Macdonald functions 



bzK {zi) h iK 1 {zi) b M , 



with argument 



= (^+p2 b 2 ) l/2 ) z _ 



[ku)b 



(3.8) 



(3.9) 



(3.10) 



while the hatted Macdonald functions are defined by K v (x) = x v K v (x). 
Substituting back into (|3.7p one obtains 



l T MN {k) = - 



mm'x 2 e lkb 



4irV "fv^ 



E 



-iT 



Kx(zi) 



<jmn - 2TjKq(zi) U(mu' 'n)+ ((2v 2 -r) 7^- - Tj K (zi) u M u N 



where 



7 2 1 

vain = {kb)u M u N ~ 2(ku)u {M b N) and V = -| = 1 - ^ + 2 ^ 2 • 



(3.11) 



(3.12) 



Note, incidentally, that as expected the part of 'zjf which corresponds to uniform motion, i.e. the term —1 in the first parenthesis, 
does not contribute to 1 13.711 . 
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Finally, omitting the terms longitudinal in w', since they give zero when contracted with any polarization tensor, 
the above expression simplifies to 11 



1 T M N(k) 



mm' x 2 e lkb 



4irV jv 3 



E 



— U 5 &MN 



2v -TJ7 T) K (zi)u M u N 



(3.13) 



As in previous analogous cases [32, |42| , the exponential fall-off of the Macdonald functions for large values of 
the argument zi leads to an effective natural cut-off N- mt on the number of interaction modes I in the sum. One 
can estimate the radius h nt of the sphere in the {l 1 } space, beyond which the modes can be neglected, by setting 



{hm/Rfb 2 



1, 



(3.14) 



from which the number of contributing interaction modes is obtained 



V 
b 1 



(3.15) 



For b <C R, which is the case of interest here, one has N- lnt ^> 1, and the mode-summation can be converted to 
integration using [53] (Z > 0) 



\ £ ^ (Vz 2 + PT 2 b 2 ) - j^^K. +d Mz) , 



(3.16) 



to end up with 



A 

jv 3 



K d/2+1 (z) 

-U o a MN 



^(2« 2 -r) 7^ - vj k d/2 {z) u M u N 



where 



(3.17) 



A: 



K D mm 
2{2n) d / 2 + 1 b d 



(3.18) 



Angular and frequency characteristics. The local radiation source (|3.17p above is expressed solely in terms of 
Macdonald functions with argument z. As has been demonstrated in analogous instances before [IHH^], the 
exponential fall-off of these functions implies an effective cut-off of 0(1) in the angle and frequency dependent 
variable z, which eventually constrains the angular and frequency distribution of the radiation. Below, it will 
be shown that a similar situation arises as well with the non-local amplitude Smn- 



B. The non-local part Smn of the source 



The stress tensor Smn is defined in (|A.5|) and calculated in coordinate representation in (|2.9p . Again, since 
its trace does not contribute to the emitted momentum, one may restrict himself to its traceless part. Write, 
as explained above, Hmn + h' MN for the leading metric perturbation and substitute in Smn, while omitting 
self-action terms, to obtain for its n-th mode (in analogy to (|2.48|l ) 



c(«) ( T fJ.\ 



-Y, 



(n-l)A,B 



ti 



(0 

NB.A 



,/(0 

" NA,B 



- h (n-l) h ,(l)AB 
~) ll AB,N n 



,M 



t MN Utl ~ 



(n-0 
MA,NB 



(n-0 
NA,MB 



(n-0 
AB.M N 



, (n-0 \ h ,(l)AB 
MN.AB 1 11 



{h <— > h'} . 



(3.19) 



and with S^ N (x^) the four-dimensional Fourier-transform of S^ > N (k ,i ). Upon substitution of the solutions 



(n) 



3.ip we have in momentum representation 



12 



? (n) 
'MN 



(k) = 



i(kb) 



4ir 2 



u M u N {ku') 2 J n {k) + 2 7 (ku')u (M J N) {k) + llJl IN {k) 



(3.20) 



11 In what follows we will omit from the radiation source all terms, which give zero when contracted with all polarization tensors, 
since they do not contribute to the emitted energy-momentum. 

12 All terms, which upon contraction with the polarization tensors give zero, have been omitted from S M ^ N (k). 
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The non- vanishing components of the scalar, vector and tensor integrals that appear in (|3.20[) . namely 

5(pu')6(ku-pu)e~ i(pb) 



I J 



if-v 2 T ) [(*-p) 2 -(4-i4-) 2 ] 



p 



(p^t), 



2 _ 1/ 

P =PmP 



(3.21) 



and 
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<5(pu')<S(fcw - P u ) e^ l(pfc) 
(p 2 -p|) [(fc-p) 2 -(fc^-^)2 

8{pu')8{ku - pu) e- i( PQ 
(p 2 -p 2 ) [(fc-p) 2 -(4-^) 2 ] 

are computed in Appendix [Bl Using (|B.8|) . (|B.15|) and (|B.17|) one obtains 14 



■Pm , 



PmPn 



(3.22) 



(3.23) 



Smn (k) 



X e Kkb) 



b b ( r 

A ^ 2 Ar 7r^T ( j/2+l(Cn) + ( 2i(ku')U( M bN) + 2i l 2v2b (M^N) ~ U M U N J #d/2«n) + 



+b 2 (tjN m Nn + ^{ku'fuMUN + 2(ku')u (M N N) j K,, , K„ ) 



1 



^^dx , (3.24) 



with C„ given in (IB.7|) . Upon substitution of Nm = — [(1 — x)(ku) + jx(ku')] um j '(^v 2 ), one ends-up with the 
fully expanded expression 



'AfiV : 



Xe 1 ^ f ( b M b N 
b 2 



o 



r yTK d/2+1 (( n ) 



i fcu') 2 ( 1 ^ I Fa f) 2 {ku) (W) 2 (1 ~ ^ (l 1 r(/, '" !J " 



(3.25) 



Like the local source (|3.17p . the non-local Smn in its effective form (|3.25j) is also confined on the brane. 
Furthermore, Smn is a linear combination of terms proportional to integrals of the form 



b 



(3.26) 



(with a = 0, 1, 2), defined in (IB. 19[) and studied in Appendix [Bf 5 . In [32| it was shown that these integrals can 
be written as a sum of two terms. One, is expressed in terms of hatted Macdonald functions of only the variable 
z and the other hatted Macdonalds of another variable z' . So, one writes in an obvious notation 



j — t [z > I n z ) 

<J(<t,t) — J ( a , T ) ~T J (<7,t) ' 



(ku)b 



, _ (ku')b 



(3.27) 



The total radiation amplitude is thus naturally written as a sum of two terms. Both expressed in terms of 
Macdonald functions, but with argument z the first and z' the second. As explained in |32| the characteristics 
of the Macdonald functions translate eventually into angular and frequency distribution properties of the emitted 
radiation. 



Note that the bulk components J™ and J" M vanish by parity only if k l T = 0. 

We passed again from summation into integration. Also, to simplify notation a bit, the index (n), labeling the number of emission 
KK-mode, is suppressed. 

Note that the coefficient A d , which appeared naturally in the definition JB . 19t of Appendix [Bj is kept here also, in order to 
avoid the introduction of more symbols. Also, note that J( a ,T) are clearly functions of k and n, which are suppressed in order to 
simplify notation a bit. 
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C. Destructive interference of the beamed high-frequency components 

As explained in [3l[ and [Hj], if the radiation source consisted of the local part Tmn alone, which is expressed 
in terms of Macdonalds with argument z, the emitted radiation would be mainly z-type, i.e. beamed $ ~ 7 _1 , 
with high characteristic frequency to ~ 7 2 /6. Now, it will be shown, exactly like in scalar radiation discussed in 
(32| . that in the sum Tmn + Smn the two leading powers of 7 in the ultra-relativistic expansions of Tmn and 
Smn exactly cancel and the total amplitude is suppressed, compared to each one separately, by two powers of 
7. This cancelation is a consequence of the equivalence principle: geodesies of the ultrarelativistic par ticles and 
of the emitted massless quanta stay close together and enhance the formation length of radiation [39( . 

We will not repeat here the derivation presented in [32j. Instead, we will just highlight the main steps of the 
derivation in the present context. 

Step 1: For z-type radiation, one has z ~ 0(1) and z' ~ 0(pj). Thus, the contribution of 1, which contains 
only Macdonalds of z' are exponentially suppressed and can be neglected. Thus, the dominant contribution to 
z-type radiation is due to . 

Step 2: The integrals appearing in Smn with a = 1,2 are suppressed compared to the one with a = 0. 
Consider for instance J(i T )- Using (|B.22lHJ.32[) and (|B.35[) from Appendix [Bl one can write for its leading terms 

J a,r) = ^ (-sm 2 #K d/2+T+1 (z) - {d + 2 ^ )p2 K d/2+T+2 {z)+ J^k d/2+T+3 (z^j + ... (3.28) 

with a, /3 and £ defined in Appendices [Al and [Bl This is of order of a~ 2 -C 1 compared to J$ T \ (given in 
Appendix |B|) . since in the region of interest a = 0(7). 

Step 3: Since the destructive interference effect has to do with the two leading terms in the local and non-local 
amplitudes, we focus on exactly these two. Hence, we neglect (ku)/~f compared to (ku 1 ), set T = 1 and v = 1 
and omit x— and x 2 — terms in the integrand of (|3.25p . Thus, the leading part of Smn reads 



S M N^Xe l(kb) 
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1 

N lK d/2+1 (( n )+ (^iu {M bN)(ku')-^UMUN)K d / 2 (C,n)+z' Z ^k d / 2 _ 1 (C,n)uMUN 



ir- 



< kb ^dx . 



(3.29) 



Step 4'- The quantities to be compared are the projections of Smn and Tmn on the polarization tensors. It 
is easy to check, using the formulae of Appendix [A] for the contractions of u M and b M with the polarization 
vectors, that at small angles < 7 _1 no extra powers of 7 are introduced through such contractions. Thus, 
comparing directly the coefficients of umUn, uim^n) and 6jv/ & jv , one ends-up with the following expression for 
the leading terms of Smn- 



1 

2 



S M N^\e l(kb) 2iu {M bN)(ku')k d/2 (C n ) + z ,2 1 k d/2 _ 1 (Cn)uMUN e-^ x dx, (3.30) 



Step 5: Using (|B.10[) appropriately simplified for the range of parameters of interest here, namely 
1 

/ dxe-^kr-iiCn) C J-k r {z)-i-^-^k T+1 {z) + ... , 

j jzz' [pfzzy 





into p.30p . one obtains 

S MN * \e l{kb) 



— K d / 2 {z)u M UN - i^^K d / 2+1 (z)aMN 



(3.31) 



plus contributions to the real and the imaginary parts suppressed by two powers of 7 compared to the terms 

kept in (j3~3T) . 

Step 6: Making the same approximations for the local part 1 Tmn i n (|3-17[) . one ends-up with 1 Tmn = ~Smn 
given in (I3.31|) above and proves the so called destructive interference of the leading terms of local and stress 
amplitudes. 



D. The total radiation amplitude 



To summarize, the total radiation source tmn was naturally split into the sum of a function of z and a 
function of z' . The two leading terms in the expansion of the function of z in powers of 7 were shown to cancel. 
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The next step is to compute the leading terms of tmn-i which survive destructive interference and use them 



to compute the emitted radiation energy to leading order. The discussion follows the steps described in 32 1 
for scalar radiation, only repeated three times because of the three relevant polarizations. That is, for each 
polarization one first distinguishes three characteristic regimes of angular and frequency distribution of the 
emitted energy, namely (a) the beamed ($ ~ I/7) high-frequency (cj ~ 7 2 /6) regime, (b) the beamed medium- 
frequency (w ~ 7/6) regime, and (c) the unbeamed ~ 1 medium- frequency regime, and then computes the 
emitted energy in each one of them. 



1. The beamed, high-frequency radiation amplitude 

In this regime z' ~ 0(7) and the contribution of Macdonald functions of z' is exponentially suppressed. 

Polarization I: Projecting (j3~TTl) and (j3~25l) on sf iN and taking into account that only the e^e^ /\/2 part 
of ef iN contributes, one obtains for x Ti = l T MN e^ and Sj = S MN ef N 



l Ti(k) 



X e i(kb) 



-iT 



Kd/2 -^(kb)[^ 1] + ((2v 2 -r) 7 ^ - r) i^ /2 (z) 7 V sin 2 



(3.32) 



and 



S 



A7 e 



(fefc) 



V2v 

( 

xK d/2 (( n ) + 



(cos 2 0cos 2 - sin 2 (f>)TK d/2+1 (( n ) + 



2iv cos d(kb) 1 1 - S - ^ T{1 5 x) 



-Tsui 2 !? 



If 2x Tx 2 \ Jl-x)f Tx\ r(l-x) 



-0 2 y v 2 v 



ipv 2 



z 2 v 4 sin 2 #K d/2 -i(Cn) \ e-^ x dx , (3.33) 



respectively. The next step is to expand in powers of 7. For that one has to use the expansions of cos$, v, T 
and ip in Appendix [5] the fact that integrals with extra x a in the integrand are suppressed by 7 _2a , and the 
expressions for the integrals jfe' > and J$ T \ given in Appendix iBl and (I3.28[) . After a somewhat tedious but 
straightforward calculation one ends-up with the leading part of n = tmn^ 1 ^' 



n ~ 



A e^l sin 2 ■& 
y/2 7^ 



d + 



d+2 



X V („\ ■ I A 1 ^ ^/2+l( Z ) , f 7' 
-•Kd/ 2 (z) + ^ + I _ 



lb 2 cos 2<A . 2 
- sin 1 



sin 1} 



K d /2+2{z) 



(3.34) 



Polarization II: Similarly, projection onto the second polarization tensor leads (with an analogous notation) 



to 



and 



, 2»Ae'( fcb > K d/2+1 (z) 
±11 = = 7 sin t7 sin 0. 

V2u z 



(3.35) 



S n = 



7A e 



i(kb) 



V2v 



r cos sin 20K d/2+1 (( n ) - 2iz sini? sm<£ ( — ( 1 - -j ) - T(l - x) J K d/2 (( n ) 



Tx 



e- l(kb)x dx . 



Following the same steps as above, one obtains for the leading ultra-relativistic term in their sum 

t\\ ~ =— 7(w-|-sin #)sm2<2> „ . 

V2 2^ 



(3.36) 



Polarization III: Repeating the same steps for the third polarization tensor, while taking into account that 
only the 6m bjy— binomial does not vanish upon projection on e 2 , one obtains: 



d t 



d+2 



Ti, 



(3.37) 



and 



5, 



in 



d + 2 



Si 



A7 e 



i(fefc) 



V2 



2 sin 2 (bK d/2+1 (Cn)e~ likb)x dx 



(3.38) 
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Combining these two and expanding, one has to leading order 

-*< /9 (*) + (d + !)__+(_ 



Till : 



d A e*( fc6 ) sin 2 19 
d + 2 7^~ 



V' 2 



d + 2 



sm 



^d/2+2^) 



Up to the common phase factor and the factor A in front, all three amplitudes are real and of 0(1/7). 



2. Wide-angle, medium-frequency radiation amplitude 



(3.39) 



In this regime it is more convenient to consider directly the behavior of the radiation source, before one 
projects onto the three polarization tensors. Also, in this regime z ~ 7, and consequently the Macdonald 
functions of z give exponentially suppressed contributions and will be neglected. Thus, the leading contribution 
is due to the z'— dependent part of the non-local source Smn, which we compute next. 

In Subsection 13 CI it was shown that for lu ~ 7 2 /6, 1? ~ 7 -1 one has J(x, T ) ~ J(o,t)/i 2 - This implied that only 
the neighborhood of x — contributes to these integrals. 

Similarly, it will be shown that in the region of interest here, namely (u ~ 7/6, 1? ~ 1), the leading contribution 
to these integrals comes from the neighborhood of a; = 1. Indeed, consider the integral with (1 — x), i.e. 
J(o, T ) — J(i,t)- Using (|B.32[) and subsequently (IB.22|) this difference takes the form 



Ji 



(0,r) 



J{l.r) 



2\ K d/2+T+1 {z') cos 2 p 

1 7 9 1 ~ A d/2+T+3 



(*') + 



(3.40) 



This is of order C(7~ 2 ) with respect to both J(o, T ) ancl J(i,t)- This implies that the integrals J(o,t) an d J(i,t) 
are almost equal, and the main contribution in them comes from the domain x = 1 — 0, and consequently, that 
integral with (1 — x) 2 is even more suppressed. 

Thus, it is natural to rearrange the terms of the integrand in (|3.25|) in powers of 1 — x as well as in powers 
of 7, taking into account that each power of 1 — x contributes a factor of C(l/7 2 ) to the integral. Taking in 
addition into account that in the region of interest here "f(ku') is of the same order as (fctt), one obtains 



ri(v , Ae W/|^ /2+l(Cn)+ 



b 2 



, , /,;>, (ku)\ , > d+l(ku')\ u M u N 
2iu {M bN)[[(ku) — ){1-X)- 



d + 2 7 2 



-u M u N K d/ 2-i(tn)\e- imx dx. 



d + 2 7 

Using (e±u) = jvsm-d, (e 2 u) — and {e\b) = 0(b) = (e 2 6), one can simplify (13.41)) and write instead 



K d/ 2(Cn) 

(3.41) 



T MN - 



7 



^M^N 2t> (/- \ i> If \ 

b 2 7 Kd/2+l{Cn) ~ U M U N K d/2 (( n ) 



1 2 

z u M u N K d/2 _ 1 (C n ) 



d + 2 



-r{kb)x dx _ (3 42 ) 



It is straightforward to check that its projection onto any of the three polarization tensors gives to leading order 
the same result as the projection of (I3.41[) . 

From (|B.22[) one obtains to leading ultra-relativistic order 



< kb >K T {( n )dx 



e — i(fcb) 
z /2 'f 2 ip 



k T+1 {z'). 



using which in ()3.42j) . one finally ends- up with 



A 



T~MN 



z' 2 7 3 ?/> 



bM ^ N 7 2 K d/2+2 (z') - u M u N k d/2+1 (z') - 



d+2 



z' 2 UMU N k d / 2 (z') 



(3.43) 



(3.44) 



3. The beamed, medium-frequency radiation amplitude 

Starting again from (|3.17[) and (|3.25|) . one can observe that for ~ j^ 1 the projection on polarization tensors 
does not introduce 7-factors, hence one can estimate the order in 7 directly from the coefficients of the tensor 
binomials umun, &m^JV and u^m^n)- I n t ne regime of interest here one has z ~ I/7, z' ~ 1, using which 
one can immediately conclude that both local and stress sources, as well as their projections on the relevant 
polarizations are of the same order, namely Tmn ~ Smn ~ t ~ 0{"f) in any dimension. 
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E. Summary 

The results of this section are summarized in Table I below, which shows the leading behavior of the corre- 
sponding amplitudes after projection on the polarization tensors. 



&\ 


U3 ~ 1/6 


LJ ~ 7/6 


cj - j 2 /b 


I- 1 


no destructive interference 


no destructive interference 
Si*] ~ T ~ ~ 7 


destructive interference: T ~ — S'M 

5M~exp(- 7 ), 

r = 0(T/ 7 2 ) ~ I/7 


1 


no destructive interference 


destructive interference: 

w T ~ exp(- 7 ) 
t = 5 = S"! 2 '! - 7- 1 


destructive interference 
T r~j S ~ r ~ exp(— 7) 



Table I. The leading behavior of the amplitudes (projected on polarizations) in the various characteristic angular- 
frequency regimes. S 1 ' 2 ' and S 1 ' 2 ' stand for the contributions of integrals and j' z respectively. The dependence of 
A on 7 is not included, i.e. all estimates are given in units A = 1. 

Notice that the entries in Table I do not depend on the dimensionality d. Also, note that this leading 
dependence of the amplitudes on 7 is valid in the whole range [0, 2ir] of the azimuthal angle tp. 



4. THE EMITTED ENERGY 



The emitted energy, obtained from (|2.43p for M = is 



E = 



4(2^-1 ^ 



J Lu D - 2 duj J dn \T LR {k)e^ R \ 

S D-2 



(4.1) 



It differs from the scalar radiation case studied in [32] by the summation over polarizations and the substitution 
of the scalar charge / — > x^m. Thus, it is expected to lead up to numerical coefficients of order 0(1) to the 
same behavior in 7 as the energyemitted in the scalar case. On the basis of the qualitative arguments given 
in the beginning of Section 4 of [32J and the subsequent numerical computations one obtains 16 for the emitted 
energy in the various angular-frequency domains the behavior summarized in Table II. 







Wfl ~ 7/6 


uj d ~ J 2 /b 


lod > 7 2 /b 


7- 1 


negligible 
(phase space) 


E ~ 7 3 , from T and S 


E~-f d+2 , fromT + S'W 


negligible radiation 


1 


negligible 
(phase space) 


E~Y +1 , frornSl 2 '] 


negligible radiation 


negligible radiation 



Table II. The relative contribution in the total emitted energy from different characteristic regions of angle and 
frequency. All estimates are given in units A = b — 1. 



A. The total radiated energy 

According to Table II the frequency and angular distribution of the dominant component of radiation depends 
on the dimensionality d. Thus, like in the scalar case [32j the study of the emitted energy has to proceed 
separately for d 2 and for d = 0,1. Furthermore, the following analysis is a triple repetition of the scalar 
case, because of the three polarizations. 



Up to a logarithmic overall factor for d = 1, discussed below. 
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d ^ 2. 

For d ^ 2 the dominant radiation is beamed with characteristic frequency around uj ~ 7 2 /&- According to 
Table I, the emission amplitude is C(7 _1 ), dominated by its real part. 

For R 3> b the KK-summation is replaced by integration and the relevant formula for the emitted energy is 
(|2.43p . One next substitutes the high-frequency amplitudes (|3.34l 13.361 and 13.39)) corresponding to the three 
polarizations and neglects their behavior at smaller frequencies. Squaring each amplitude, one recognizes the 
appearance of products of K d / 2 { z ), Kd/2+1/ z 2 an d Kd/2+2/ z 2 with themselves, and also the presence of the 
factor uj d+2 from phase space. Thus, one starts with 



dE 
~dh 



4(2tt) 



D-l Z-^/ 



D-2 



du \T L R{k)e v 



LR I 



(4.2) 



and upon integration over the angles except $ one obtains 



dE 



(^mm')V +1 sin^tf ^ ^(d) n (d), 



8 (27r) 2d + 5 b 3d + 3 ip d + 3 



a,b=0 



where 



(d) 



J K d /2Uz)K d /2 +b (z)z d+2 ^+ s ^- 1 ' 



dz . 



(4.3) 



(4.4) 



(d+l) 3 sin 2 tf 
(d + 2) 3 7V 2 
d + l 



D<?0?) = (d + 2)£>M 
Z?ff0?) = (d + 2) 2 £>$(tf) 



(d + 2) 2 



, , , , sin 2 $ 



d+l f 2 T^p 2 , 3sin 2 i? 



d + 2 V sin 2 tf 4 7V 2 



d + 2 



+ 7 2 sirTtf . (4.5) 



(all of O(l) for ~ I/7), while use was made of (|C.5D : 

J dfl d+ i = Q d+ i , J cos 2 (j}dQ d +i = Tjf^d+i , J cos 4 <j}dQ d +i = ~^d+i 

Upon integration over $ using (|C.2I) one ends up with 



E = 



I (27r) 2d + 5 r(d/2 + 1) 6 M + 3 



a, 6=0 



where the constants D a d J are given by 



' d+6 \ p /d+4^ 
(d + 2) 3 T(d + 5)' 



n (d) g d +5 (^+i) 3 rm r m 



DW = (d + 2)D, 



mi 



D, 



(d) _ _ d+3 



(d + i)(2d + 3)r(^)r(^ 2 ) 



(d + 2) 2 



r(d + 4) 



D$ = (d + 2) 2 I? ( 



2 n M 
00 



(4.6) 



(4.7) 



(4.8) 



D[f = (d + 2)D, (d) 



112 



£)W = 2 rf - 2 (5d + 6)(3d 2 + 15d + 20) 



r(d + 4) 



Finally, after the summation in (|4.7I) one obtains for the total emitted energy in the lab frame 

(x 3 D mm') 2 d +2 

g ~ Cd b 3d+3 7 

with C 2 - 6.23 x 1CT 4 , C 3 = 2.61 x 1(T 4 , C* 4 = 1.98 x 1(T 4 , C 5 = 2.16 x 1(T 4 , C* 6 = 3.10 x 10~ 4 . 



(4.9) 



d=l. 

The emitted energy flux as well as its frequency and angular distributions in d = 0, 1 were computed numer- 
ically. 
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<H ■ 1 ■ 1 1 1 ■ 1 > 1 1 1 1 1 1 1 r 

0123456 2468 

cob/y Q-J 



(a) (b) 
Figure 3: Frequency (a) and angular (b) distribution of the total emitted energy for d = and 7 = 10 3 . 

According to Tables I and II in this case most of the radiation is beamed d < 1/7, but, as it was described 
in [32( 1 and immediately follows from the destructive interference, in five dimensions the frequency distribution 
falls as 1/uj in the region from 0(7/6) to C(7 2 /6), hence the total emitted energy grows as 7 3 ln7, while the 
dominant contribution comes from this entire region. From the known expressions for the amplitudes in integral 
form, we compute the integral numerically for several values of 7 and conclude: 

E~C 1 ( ^T l/)2 7 3 hi7, d = 1.66 x 10" 4 . (4.10) 



d = 0. 

Again, according to Tables I and II in this case most of the radiation is beamed within $ < I/7 and with 
frequencies around uj ~ 7/6, with the contribution from higher frequencies decaying as 1/uj 2 . In this domain 
the local and non-local parts of the amplitude are equally important. Thus, one has to add (|3.17p and (13.251) 
to compute the total radiation amplitude, project onto the two polarizations e±, square each one of those, add 
them up and multiply by the appropriate phase space factor. After the angular integration one obtains the 
frequency (or z') distribution of the emitted energy shown in Figure [3^, while the angular distribution of the 
total emitted energy (which carries the anisotropy) is given in Figure 

Integrating the frequency distribution numerically over w from to +00 one obtains for the total emitted 
energy 

£^o W 7' )2 7 3 , C = 5.7xl0- 4 . (4.11) 



B. The low frequency part of the spectrum 

A few comments, related to the low frequency part of the classical radiation spectrum, are in order here. 

(a) For d — the distribution dE/dui of the emitted energy has non- vanishing finite zero- frequency limit. 

Indeed, comparing local (|3.17l) and stress (|3.25[) sources one concludes that for ui — > +0 the dominant part of 
the amplitude is given by the imaginary part of the local piece 17 (omit the phase factor e 4 ^^), whose behavior 
is 

XTK^z) 1 

TMN\k) ^ —t o ~ — ■ (4-12) 



In the diagrammatic quantum mechanical treatment of soft gravitational radiation there is a corresponding statement, namely 
that the soft radiation emitted from internal lines is negligible. Only gravitons attached on the external lines contribute in the 
zero-frequency limit. 
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Upon contraction with the polarizations e±, substitution of the finite zero-limit Ki(0) 
Macdonald function, and setting v = 1, F = 1, one obtains 



iX 7 sin i9 
~2 tub 



1 



7 2- 2 



1 \±2i- 



1 of the hatted 



(4.13) 



Substitute this into (|4.2[) and integrate over the azimuthal angle <f), to obtain 



dE 
duj 



16(2tt) 2 6 2 



dtf sin 3 tf 1 + — 



1 



i/> 2 ' 7 4 v> 4 



(4.14) 



plus subleading powers of 7. Note that in the low frequency limit all angles (0 ^ ■& ^ ir) contribute to dE/duj. 
Finally, using V$ -4/3 f03j) . U 4 3 ~ 4 7 4 /3 JUSJ, U 2 3 ~ 4 (In 27 - 1) (|04|) . and substituting the expression 



for A, one ends-up with 



dE 
dcu 



j=0 



(xrfrnm 
4(2tt) 4 6 2 



7^(ln2 7 -- 



(4.15) 



which for yc\ = m = ml = b=l and 7 = 10 3 agrees with the value obtained numerically in Figure [3ji 18 . 

(b) One can multiply the above result for dE/dw\ u —o by the range w max ~ 7/6 to obtain a rough estimate of 
the total emitted energy. Its 7 dependence differs from f|4. by a logarithmic factor. 

However, this simple method to obtain an order of magnitude estimate cannot be used in d > 0, because 
according to Table II the dominant contribution to the emitted energy comes from the high frequency regime, 
which does not have overlap with the u) — > domain. 

(c) For d > 0, one cannot replace at small frequencies the KK-mode summation with integration. One has 
instead to keep separately the three-dimensional coordinates w, 0, if and the extra-dimensional discrete momenta 
kf- In this case 



(ku)b 
7« 



kj, — vjv cos 9 



and 



(ku')b 



= — \lm 2 + k% , 



while (kb) = — njsinflcos^. 

The real part inside the square bracket of (|3.13|) is regular or blows-up logarithmically at vj = 0. Thus, 
it vanishes when multiplied by the volume measure. Consequently, the dominant contribution to the emitted 
energy is due to the imaginary part, i.e. (ignoring the irrelevant phase factor) 



TMJv(fc) ^ 1 T M N{k) 



■D 



Kx( Zl ) 



-&MN ■ 



(4.16) 



Taking into account that K\(z{) is regular at all frequencies, the low- frequency behavior is determined by the 
factor <tmn/z 2 — lo/z 2 . For k\ > this limit is zero and does not contribute. For kx — we have w = oj 
and thereby ctmn/z 2 ~ 1/w as in four dimensions. Thus, as expected, only the zeroth emission mode n = 
contributes in the summation over n of (|2.53p in the uj — > limit. 
Now, for kr = z — usb4>(0)/v with ip{6) = 1 — v cos9 and 



TMN(k) 



mm! x 2 



— 2^ 9,.9.,.9//n a MN , 



Aw-yV ^ uj 2 b 2 %b 2 {9) 



Pt 



V 
R 



(4.17) 



For b <C R, which is the case of interest here, the summation over the interaction-KK modes can be restricted 
up to modes with \pr\b — 1. The result, which can also be obtained by the usual substitution of the summation 
with integration, is 



... A 2 d / 2 T(d/2 + 1) 

T MN\k) — .9 L 9_,.9/m °~MN 



7 L0 2 b 2 lp 2 {9) 

and subsequently, using (|2.53p . one concludes that 

2iJ-2r2w/oj.iVJ! .,.,.,<•*■! 

u=0 



dE 
duj 



2 T 2 (d/2 + l)(>c i D mm')' , 
(2n) d + i Vb 2d + 2 7 



In 27 



1 



6 V d+2 



(4.18) 



(4.19) 



Note that this result differs from the one obtained in [lql . where it was claimed that dE/dw blows-up logarithmically as la ■ 
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The volume V of the internal space survives in the denominator, because, as explained above, only the n = 
mode of radiation contributes. Of course, for d = f|4. 19[) coincides with (|4.15j) . 

(d) In M £)>4 the frequency spectrum vanishes for uj — > 0. There are many ways one can convince himself 
about this. First, it is obtained by taking the limit V — > oo of (|4.19l) . Alternatively, observe that in higher 
dimensional Minkowski the relative contribution on the various radiation modes satisfies En=n /En 2 >n — ^ 0- 
Taking in addition into account that non-zero modes n do not contribute in the limit uj — > 0, one concludes 
that dE/duj \ ^=0 vanishes in higher dimensions. Finally, the same conclusion is reached by taking the behavior 
r ~ square and multiply by the measure ui d+2 . Again, for d > the limit of the integrand for to — > is 
zero. 



C. Quantum constraints 

The above results were obtained in the context of classical theory. However, quantum mechanics restricts 
the region of validity of a classical computation. Thus, to ensure the reliability of the above conclusions, apart 
from the condition b <C R, which allows one to replace the summations over KK modes with integrations, and 
which will be assumed in the sequel, one has to justify (a) the use of trajectories for the colliding particles, 
(b) the expansion around flat space for the gravitational field, and (c) the use of classical field description of 
the radiation. From the analysis of the elastic ultrarelativistic scattering problem [42| one concludes that the 
classical computation around flat space-time of the cross-section is identical to the eikonal approximation of the 
quantum answer, as long as the conditions 

y/i = 2m 7cm » AU , b > 6 7 = r s7c ^ d+1 ) , b < b c = r s (r s /\ B ) 2/d (4.20) 

are satisfied 19 . For the radiation problem at hand, the weak particle- recoil condition due to the emission of 
gravitons with characteristic frequency u> is satisfied if the momenta of the emitted gravitons are much smaller 
than the momentum transfer of the elastic collision. However, experience from analogous computations of the 
total energy of synchrotron radiation shows that this condition can be relaxed and replaced, instead, by the 
weaker uj <C £q, where £q — fni is the total energy of the colliding particles. When the emitted energy E is of 
order £q, this condition also guarantees a large number of emitted quanta, and justifies further the description 
of radiation with a classical field 20 . 

1. The lj ~ J 2 /b and uj ~ 1/b regimes 

• Radiation with characteristic frequency uj ~ 7 2 /b does not satisfy the above constraints, and consequently 
may not be reliably treated classically in the presence of extra dimensions. Indeed, for uj ~ 7 2 /& the condition 
uj < my leads to b > 7/m > b ci which contradicts the requirement b < b c . Thus, in the sequel we shall ignore 
the radiation with characteristic frequency 7 2 /6, even though according to Table II it is classically dominant 
for d > 2. 

• On the contrary, as explained in [32J and shown in Table II the radiation emitted in the soft frequency 
regime uj < 1/6, satisfies all the above conditions, but its energy is negligible by phase space in all dimensions, 
and will not be discussed further. 



2. The uu ~ 7/6 regime 

According to Table II the dominant radiation in this regime is beamed for d = 0, 1, almost isotropic for d ^ 3, 
while for d = 2 the two components are of comparable importance. 
For uj ~ 7/6 the quantum condition w < 7717 reduces to mb > 1. 
Thus, in terms of the rescaled dimensionless parameters 

A^^, P = bM* (4.21) 



For colliding particles with equal masses 7^ m = (1 + t)/2 ~ y/2 

In what follows and in order to simplify notation we replace the strong inequalities by simple ones. 
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a reasonable set of requirements for the validity of the above classical results in this intermediate frequency 
regime is 21 : 

£2\ i/(<m) 



£ > 1 , /i « 1 , > ( — J , /3 < £ 2/rf , > 1 , (4.22) 

which are equivalent to 

1 /£2\V(d+l) 

£ > 1 , m<1, - < ( — I </?< f 2/d , (4.23) 

M \ M / 

and can easily be satisfied in all dimensions 22 . 

Thus, we shall next compute the energy emitted in the uj ~ 7/6 characteristic frequency regime 23 . 



D. Radiation emission with characteristic frequency ui ~ 7/6 

As explained above, the estimates in Table II were obtained under the assumption that the frequency inte- 
grations could be restricted to the appropriate regime corresponding to each entry, i.e. from a fraction of the 
corresponding characteristic frequency to a few times that frequency. The quantum conditions discussed above 
restrict lu to lu ^ u) max = mj, which translates to z' mayi = mb > 1. 



1. d > 3 

Thus, for d ^ 3, according to Tables I and II the leading contribution to the emitted energy is obtained 
from (|3.44[) integrated over phase space with z' up to z' ma ^. But, the fact that the integrand is exponentially 
suppressed beyond z' ~ 1 , allows one with negligible error (in the overall coefficient) to extend the z' integration 
to +00. As far as the lower limit is concerned it is a priori 1/6, since Q3.44p is valid for u > 0(b^ r ). However, 
since the integrand in (|4.4p is regular at z' = one can shift the lower limit of integration to 0, the relative 
error being of order I/7 2 . This reproduces the power of 7 estimated in Table II. The purpose of the rest of this 
subsection is essentially to obtain the overall coefficient. 

Projecting (|3.44[) onto e± and em one obtains 



X 



2z' 2 jip 



(cos 2 4> cos 2 t9 - sin 2 </> ± i sin2cj)cosi!})k d /2 + 2(z r ) - sin 2 i?( K d / 2+ i{z') + 



d + 2 



z' 2 k d/2 { z ') 



A 



2(d + 2) z^-yip 



(1 - cos 2 sm 2 d)K d/2+2 (z') - sm 2 d[K d/2+1 (z') + 



d + 2 



z' 2 k d/2 (z>) 



(4.24) 



One then according to (|2.43[) has to square these, multiply with the phase space factor and integrate over angles 
and frequencies. 

The result of the integration over frequencies takes the form 



(x^-mm 1 ) 7 



n2„,d+l 



dE 

dfl d+2 ~ 16 (27r) 2d + 5 b 3d + 3 V* 2 



(4.25) 



where C a ^ are given by (14.41) and 



D$(d,<l>) 
Df 2 \^^) 



1 



(d + 2) 3 
sin 2, i? 



d + 2 



cos 



-(d + 2)D$(d, 



d + 2 



(cos 2 sin 2 i? - 1) 



4?(0,0) = (d+2)£W(0,0 



D^{d,4>) = {d + 2fD^{§, 



d+l 
d + 2 



(1 - cos 2 </> sin 2 i?) : 



(4.26) 



21 The condition fi < 1 is necessary to justify the use of point-particle, instead of black-hole description of the colliding particles. 

22 For d = there is no upper bound on the allowed values of /3. 

23 Incidentally, this is the frequency regime in which one would a priori expect to have dominant radiation: Due to the Lorentz 
squeezing of the field lines along the direction of relative motion of the colliding particles, the characteristic duration of the 
collision in the lab frame becomes r ~ 6/7. Thus, the characteristic frequency of the emitted radiation should be uj ~ 1/r ~ 7/b. 



20 



The values of c£i are to be computed with help of (|C.6[) . 

The integration of (|4.26[ [4725]) over the angles except -d (i.e. over the sphere S d+1 ) is easily performed using 
formulae (|4.6p , while the final integration over ■d is done with the help of (|CJ.3[) 24 . Making use also of (|4.4p , one 
finally obtains for the emitted energy 

- (x^mm') 2 d+1 
E - Cd b 3d+3 7 . I 4 - 27 ) 

with <5 3 = 1.08 x 1(T 4 , C 4 = 4.77 x 1(T 5 , C 5 = 4.03 x 10~ 5 and C 6 = 4.98 x 10~ 5 . 

2. d= 

As explained above, the z' integration extends a priori up to z' max ~ m& > 1. It was checked numerically that 
(a) the estimate of the cubic power of 7 given in Table II is correct, and (b) that the coefficient is insensitive to 
the upper limit of integration, varying by a factor of two as one changes the upper limit of the z' integration 
between 1 and 00. 

Thus, (|4.1ip as a reliable estimate. 

3. d=l,2 

In this case the estimation of the contribution to the emitted energy from the oj ~ 7/6 part of the spectrum 
was done numerically. Specifically, one started with fl3.17[) and (13.25)) and integrated numerically over 1} from 
to 7r and over z' from to something of order 1. For a reliable estimate it is enough to integrate up to z' = 1. 

The result for d = 2 is 

ff^O.04 ^g^V In 7 (4.28) 
b J 

and for d = 1 

E ~ 0.9 x 10- 4 M^l! 7 3 . (4.29) 

E. Extreme radiation efficiency for d ^ 3 and uj ~ 7/6 

Consider specifically the case m = to', divide the emitted energy by the total energy £0 = m + m-f of 
the colliding particles and use the definition of r$ given in the Introduction, to write the radiation efficiency 
e = E/£o in the form 25 

e~(-£) 7^-7^7^, d = 0,l (4.30) 

and 



6/ M /?3(d+i) 



e~(fj 7^- ^3(^+1)^-3 . ( 4 - 31 ) 

It should be pointed out that these expressions differ from the generic formula e ~ (?'s/fe) 3< - d+1 - ) 7cm +1 given 
erroneously in [33|. The difference is due to an error in 33] related to the extent of destructive interference. 
Nevertheless, the qualitative conclusions presented in [33[ are still valid. Incidentally, the powers of 7 in the 
emitted energy differ from the scalar radiation studied in [32[ only by the substitution / — > X£>m. 



7T 7T 

24 At small angles expressions 114.241 1 are not valid. But as before, if we substitute integral J dd by J d9, the relative error is 

0(1/7) 

0(1/7) since the integrand doesn't blow up at $ = + 0(1/7). Thus we may apply JC 3l t ■ 

25 Since it is not our intention to engage in numerology, we ignore in these formulae inessential d-dependent coefficients, even though 
they are roughly of O(10 d ), or the factor ln7, which is present for d = 2. 
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In the parameter range denned in (|4.23l) the radiation efficiency is always much smaller than one for d = 0, 1, 2. 
Indeed, consider the case d = 0. The maximum e is obtained for the minimum allowed value (3 ~ (£ 2 //i) 1 /' <i+1 ) 
of the impact parameter and is equal to 

Cmax(d = 0) ~ (|) « 1. (4.32) 

Thus, in d = a small fraction of the available energy is emitted in gravitational radiation and vanishes in the 
massless limit. However, one can convince her/himself that for d 3 the efficiency can take values arbitrarily 
close to one or even become much greater than one for a wide range of parameters obeying the classicality 
conditions (|4.23l) stated above. For instance, in d — 4 the minimum value of e is obtained for (3 ~ £ 2 l d and is 

e min {d = 4) ~ ^ > 1 . (4.33) 
A 1 

Ditto for dimensions d > 4. This is a priori a puzzling result. It leads for d 3 to unphysical, i.e. greater 
than one, values of the radiation efficiency, and actually blows- up in the massless limit (/i — > 0). In relation 
to this we would like to offer the following comments: (a) The massless limit is a very singular limit in our 
approach. The characteristic frequencies and angles all become singular in that limit. In particular, since in 
that limit 6 c /6 7 ~ (£ 2 / d fJl ) 1 /( d + 1 ) -j. 0, it is impossible to satisfy P~2"51) . even if one ignores the "quantum" 
constraint /3/x > 1. Thus, the massless limit lies outside the domain of validity of our approximation, (b) 
Weinberg has shown [4l[ that in d = there is no infrared divergence in the soft graviton emission rate from 
massless colliding particles. This is in agreement with our results not only for the soft component in d = 0, 
but also for soft as well as high frequency emission in d = 0, 1, 2. (c) The explicit computation presented here 
is rather involved, but, as explained, the powers of 7 in the emitted energy are determined on general grounds, 
based on properties of Macdonald functions, number of KK interaction and emission modes and dimensionality 
of phase space, (d) One may suspect that some classicality condition is missed, which might render the above 
formulae non-applicable. But the constraints used in (I4.23[) are quite straightforward and based on well-known 
physics. So, (e) it may well be the case, that the above results are reliable and the interpretation of their 
puzzling features is that in ultra-planckian particle collisions there is a lot of energy emitted in gravitational 
radiation. This will lead to strong radiation damping, which should be taken into account in the study of such 
ultra-planckian scattering processes. 



5. CONCLUSIONS 



A detailed study was presented of classical gravitational radiation emitted in ultra-relativistic collisions of 
massive point-particles interacting gravitationally. The space-time was assumed to have an arbitrary number 
of toroidal or non-compact extra dimensions and the post-linear approximation scheme of General Relativity 
was employed for the computation. The angular and frequency distributions of radiation, as well as the total 
emitted energy were studied in detail to leading ultra-relativistic order. 

Three characteristic frequency regimes (1/6, 7/6 and 7 2 /6) of the emitted radiation were identified and the 
characteristics of the dominant contribution was determined in various dimensions. 

In particular, in any number of dimensions the soft component of radiation is mainly due to the scattered 
particles, with negligible contribution coming from the cubic graviton interaction term 26 . However, in contrast 
to the four-dimensional case, in any number of extra dimensions d > the frequency spectrum of the emitted 
radiation vanishes as uj — > and the total emitted energy in soft gravitons is negligible. 

Also, contrary to what happens with the soft radiation emission, the cubic graviton interaction and the 
scattered particles themselves are equally important as sources of radiation with high frequency. In fact it was 
shown that in any dimension they lead to partial cancellation (destructive interference) of the total beamed 
radiation amplitude in the high frequency domain, as a result of which the emitted energy in the 7 2 /6 frequency 
regime is reduced by two powers of the Lorentz factor 7 in the Lab frame. 

The relevance of the classical analysis to the full quantum radiation problem was also discussed. The classi- 
cality conditions, necessary for the classical treatment to be a good approximation to the full quantum problem 
were derived and the radiation efficiency e, i.e. the fraction of the initial energy which is emitted in gravita- 
tional radiation, was computed for parameter values inside the region of validity of our classical computation. 
Although for d = 0, 1, 2, e is smaller than one in the whole allowed range of parameters, it can be arbitrarily 
close to one or even exceed unity in d ^ 3. One possible interpretation of this "unphysical" result is that indeed 



This is a well known fact, verified easily also in the context of Feynman diagram infrared graviton summation. 
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in such ultra-planckian collisions there is a lot of gravitational radiation, which will lead to strong damping, 
which should be included in a reliable treatment of the scattering process. 

However, the implementation of this interpretation, the proper treatment of the massless limit, the extension 
of the region of validity to impact parameters outside the range (|4.23[) and the comparison with quantum results 
based, for example, on string theory, are currently under investigation and, hopefully, will be the subject of a 
future publication. 
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Appendix A: Definitions and variations 



1. Perturbation theory variation over gravitational constant in flat background 



We use the notation of Weinberg [37[ • 



9MN — n MN 
-.(l)JV 



- MN 



h, M /2 



The definition of Ricmann and Ricci tensors: 



ft NLR — 1 



NR,L 



NL.R 



r (l)-R _ ( h R 

r (l)M NR _ n 
1 NR 'I — u ■ 



h R 

a N.M 



^MN )/2 



1 SL l NR 



rJif r s 

1 SR l NL ' 



RMN = R L MLN 



(A.l) 
(A.2) 



Thus 



R 



(i) 

MN 



G 



i(»S, 

-Dh + h 



MR + h M,NR - D h-MN - fr.MJv) 



-Oh 



MN 



MN 



MN D MN > 

- K llMN 



(1) 
MN 



^(-OlpMN - VMN(,£ L + S.M.N + 6v,Af) 



--j^^MN 



£,M = d N i>MN 



(A.3) 



The second-order variations: 



i(2)Af 
LR 

,(2)£ 
ML ~ 



1 



h MN {h LN R -f h,RN,L — ft-L_R,Jv) 



1 h LR h 



r (2)M NR _ n 
1 NR 'I — u 



R 



MAT — — - ( h M B (flNB,A — hjsrA,B) — — h AB m Hab,N + h AB (HmA,NB + ^NA.M B — h A B,MN — ^MJV.As) 



R (2) = (g LR R L R) {2) = -h LR Rf R + V LR Rf R = \ ( 2h AB o hAB h AB , L h AL > B + h A B, L h AB > L 



G 



(2) 



R 



(2) 



MN ~ MN 

As a result: 



;VMN 



RW 



LR 
1 



i>M N 



RW 



(A.4) 



S 



MN 



--G {2) 



MN 



h AB (llMA,NB + h,NA,MB — h A B,MN — ^mN.Ab) ~ 7,h AE 1 M h-AB,N 



1, 

-h M N a h+ 



~^M E '{h>NB,A 



1 



lNA,B) 



,AB nU ._ u ._^AL,B + * h h AB,L 



r]MN 2h A "n h AB - h AB ,Lh J 



(A.5) 
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2. Useful kinematical formulae 



The angles in the formulae below are defined in FigH 

n l 

u' i =7(l,0,0,«), u' M = (1,0,0,0) , tp= l-ucos<9coso: = l-«cosi?, W = — 

R 

i, i zd . (ku')b ujb (ku)b uib . , , 29.22 ; l % 

vo = k , cosa=— , z=- — = — , z = - — — = — i/j = z^ip, zf = z z + bp T , Pt • = tt 

cj 'yu 'yv 'yv v ' R 

£ 2 = 2jzz' - z 2 - z' 2 = w 2 b 2 sin 2 9 + b 2 k 2 T = (w6 sin ?9) 2 = ( 7 uz' sin i9) 2 

— (fc&) = £ cos = 7z'w sin i? cos <j> = ^z'v cos a sin 9 cos </j = vub sin cos </? 

„ , 2 TO 2 & 2 COS^(l - V COS'S) 2/2,/, m , A 

P = 722' - z l = ^-5 = 7V ^(1 - ip) . (A. 6) 

v cos 2 a 

Next we compute non- vanishing products of polarization vectors: Using (j2.16[) . (|2.5ip . (|2.52[) and these 
definitions one finds the following scalar products 



(ku) = 7 w 2 + k% — vuvcos9j , (ku r ) — \j zu 2 + k\ , {uu') = 7 . (A.7) 

To calculate the scalar product {be2) one observes that e xz ° T3 - D ~ 2 b x u z u' k r e3...eD-2 = —bjvk±, where k± is 
the length of the projection of k M onto the subspace orthogonal to t, x, z, since k is transverse to all e Q 's and 
the form e T3 ' D ~ 2 k T e3...eD-2, remaining after factoring out the t,x,z subspace, represents the volume of the 
rectangle. Thus 



k± = y zu 2 sin 9 sin ip + fc^ = ui sin <j) , (A.£ 

and one finds 



, ^ 2 sin 2 0sin 2 o9 + /4\ 1/2 , /cos 2 a sin 2 9 sin 2 u> + sin 2 a \ 1/2 , . , . . „. 

(M = ~M 2 ■ 2 n , 2 =-M 5 .2.-2 =-6sm0. (A.9) 



?Z7 2 sin 2 9 + k 2 ^ J \ cos 2 a sin 2 



sin a 



Evaluation of the scalar product (e\b) is straightforward and leads to 

(be!) = - J , ro2sin2gco ^ = - * cob* (l - --^U = ^ ctg * . (A.10) 

2 J(za 2 sin 2 9 + k 2 T )(za 2 + k 2 T ) v V 7(W); (**) 



Finally, the remaining non-vanishing scalar product is 

1 /o 

tu 2 sin 2 9 + k 2 



t \ vu sm a -f re T \ / " 7-5— ; 2 ■ n /1 -1 1 \ 

(eiu)=7U 5 = 71' v cos a sm + sin a = 7usm?7. (A. 11) 

ZU Z + fcf 



Appendix B: Momentum integrals 
1. Basic scalar integral 

The local integrals for massless and massive modes are calculated in [3l| . The main integral for stress (also 
for massive modes) is calculated in [32] ]. 

^4x>^>- ^°/ \-S h ^ - (R1) 



cd = k° = \/k 2 + fc^ 2 and k is a 3— dimensional vector lying on the 3— brane, while kj* 1 = and pt* = l l /R 

with integers are d— dimensional discrete vectors, corresponding to the emission and interaction 
modes, respectively. 

Using Feynman parametrization J nl takes the form: 

jnl = f dx e - i{kb)x f d i /_M) + (ku')x]S[(pu) - (1 - x)(ku)} e-'W _ (B 2) 







[p 2 - (k T x-p T ) 2 
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Integrating over p° and splitting p into the longitudinal pu (with respect to u) and transversal pj_ parts, 
integrate over pu. The corresponding covariant splitting of 6 M on temporal, longitudinal and transversal parts 
reads 

V = ( b U ')u>»- MlM (7U * _ „M) + £ 

and denote scalar impact parameter b as 

Then, introducing in pj_ the spherical coordinates e? 2 pj_ = |p_i_|df2id|p_i_| and integrating first over the angles 
and then over |p_i_|, one obtains 

1 

jnl = ^. f dxe -KNb)K_ l{Cnl ), (B.4) 
iv J 



with 

Ql^x) = z' 2 x 2 + 2jzz'x(l - x) + z 2 (l - x) 2 + b 2 (k T x - p T f 

Nfj, = [ x (ku')u'^ - (1 - x){ku)u^\ - —2 [a;(feti')tt M - (1 - x)(fcw)w^] + xk^ . (B.5) 

Sum up over Z's with rule (|3.16[) to get 

l 

r irh 2 ~ d 

J n (k) = A d J dxe-^V K d/2 ^(C n ); A d = ^ )d/2 • (B.6) 
o 

with 

C 2 (x) = z' 2 x 2 + 2 1 zz'x(l-x)+z 2 (l~x) 2 ; Cn(0)=2, ( n (l) = z' . (B.7) 
Finally, fix (im') = (bu) — to get 

l 

J"(fc) = A d y dze-^X d /2-i(Cn) • (B.8) 



Thus 

J n (fc) = J [2] + J [z ' ] . (B.9) 
In the frequency-angular region g ~ wfr? > 1 jW and jM are given by the series over small 1/g : 

jW = ^2 (^«J/a(«) ~ i(A6)^«j/a + i(2) - (d K d/2+1 (z) + ^±k d/2+2 {z)\ + R z (B.10) 



. COStf 



i ^2j^/2+l( Z ') ) +^', (B.ll) 



where 



z 2 = 



y £ 2 sin$ 

Thus J n (k) can be split in two parts with drastically different spectral-angular behavior each. 



2. Vectorial and tensorial integrals 



Vectorial integral is denned by 



= d* P - 



5(pu')S(ku - pu)e- 1 ^ 
{p 2_ p 2 T)[[k _ p) 2_ {kir _ plr y 



-PM 



Thus 



dJ 

W 1 



Substituting J n {k) in the form (|B.6I) and differentiating, one gets taking ()B.3|) into account 



J M = K d / dxe-^ 



NMK d/2 _ x ^ n )+i^-K d/2 {Q n ) 



with 



= b M - 



(bu)u M + (bu')u' M - l[(bu)v! M + [bu')u M ] 

9 9 



N M = (N^xkr 1 ), 



where it is convenient to use the following properties of hatted Macdonalds: 



In special frame 



J M = A d / dxe~^ x 



N M K d/ 2-l(Cn)+t^-K d /2(Cn) 



Some useful products 

N -u = (ku) N ■ u' = J" • u' = 

J n ■ u — (ku)J n N ■ b — x(kb) 

Tensorial integral is defined 

J MN (fc) = I £ jj^ (fc) , ^ (fc) = / d 4 

z ^ 

Thus 



V ■ k = -r[x (z 2 + z' 2 - 2zz'-f) + (jzz' - 
1 
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V 2 = ± [x 2 (z 2 + z' 2 - 2zz' 1 ) --" 



5{pu')S(ku - pu)e~ l{ P^ 
(p 2 -p 2 T )[{k-p) 2 ~(V T ~p^y 



■ PmPn ■ 



Substituting Jjjt in the form (|B.13[) and differentiating, we have 



J MN — ^2 ' 



6m6 



mon 



b Nm Nn K d / 2 - 1 ( C" ) + 2iJV (M 6 JV) if d/2 (C„) - ^^ d /2+l(Cn)- 

1 



?7mat H — 5—^ 



UmWjv + "m u W - 2 7 U (M«W) ) K d/2((n 



Some products: 

Jmn u = \ ku )JM Jmn u = Jmnu u = (ku) J . 

Here we see that integrals J^, Jmn represent the superposition of some scalar integrals of the type 

1 





1 
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(with a = 0, 1, 2): for instance 
1 



J M — 7 
VO 



^ [z'u' M + zu M ] - [z'u M + zu' M ] + vbk M ^j J(i,-i) + z |«'m - ~ u Mj J (o,-i) + i-y-J(o,o) 



(B.20) 



Two integrals J(o,-i) an d J(o,o) (i- e - with a = 0) are of the type of the basic scalar integral and thereby known, 
while the derived one Jn t _u is new. The computation of integrals with er = 1,2 represents the goal of next 
subsection. 

3. Derived integrals for stresses 
Now consider ()B. 19|) with a = 1, r = — 1 : 

l 

J(i.-i) = A d / se-^^/a.xCGOdB, (B.21) 



o 



If we'd know J(k) exactly (IB.21|) can be calculated by differentiation: 

i dJ(o,-i) 



£ dcos()> 



Our strategy consists in the reduction of J( J i T ) int o the superposition of J(o,t') (with some t's), for which 
we already have calculated expressions (IB. 101 IB. 111) . Note that for higher r we keep and shift index of 
Macdonalds in (IB. 101 IB.lij) : thereby (|B.10I IB.lij) may be generalized into 

M Ad (av n r^MK- (<* + 2t + 1)/3 r> /3sin 2 r> \ 

J (o't-i) = ^|2 [P K d/2+Az) -t{kb)K d/2+T+1 (z) — K d/2+T+1 (z)+ fl2 K d/2+T+2 (z) I 

r 'i A^e~ i( - kbS> r o i 

Jg,r-i)~ 2.2 (^ 2 -/3)^/ 2+T (z')-^cos0^ £i/2+T+1 (z') . (B.22) 
a s L i 

Representing a; = [x — /3/£ 2 ] + f3/£ 2 , we have 



l 



•7(1,-1) = ^(o,-i) + y / re- < W a ^ lI/2 _ 1 (C„)dx = |j ( o,-i) + ^J, (B.23) 



o 

with notation 

l 

J = 

o 



J dr^e-W'Vi^)' (B.24) 

Representing [43j,|44j 



i^-i ( vV - r» ) = r^a* / ^ + ^ ^( av V+T)dy , (B.25) 



o 

we have 

oo 



J "= / d V ^l^n ^^+l) / ^e-'f^r 1 -"/,^). (B.26) 



Fix n = 1/2, then i/ = (d + l)/2 = (D — 3)/2 > -1, and 



/2 

r 1/2 /i/2 (ry)y 3/2 = y ^= sinh(ry) . (B.27) 
V 71 " 
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Taking the internal integral over x 
1 

fae^ M '>*wah(Z[x-0/Z 2 ]y) = ^ (-l) J ' _1 e- 



3=0,1 



i(kb)j ^cosh(gd J y) + i(kb) sinh^dj-y) 
ev 2 + (kb) 2 



with 



di = 3 - PI? = hi - Pit J = 0, 1 



we arrive at 



oo , 

j= Xl a (d+i)/2 /tf-i -i(kb)j f , ^(rf+i)/2(a\/y 2 + 1) y 2 £cosh£,d J y + iy(kb) smh^y 

3— o 

Substituting = — £ cos (|A.6|) to get 



oo , 

} e J y (y 2 + iy' 2 



y 2 + cos 2 (, 



In the integrand numerator add and subtract cos 2 </> to y^ 



J 



72 a ( rf+1 )/ 2 



^(-1)^-^ dy 



K( d+ i)/ 2 (uy/y 2 + 1) 

(y2 + l)(d+l)/4 



cosh £djy 



cos 2 4> cosh ^(ij-y + *y cos <f> sinh 



y 2 + cos 2 <; 



3=0,1 I 

and then, integrating cosh^djj/ in a bracket with help of (|B.25I) and comparing the remainder with 

oo 

2 1 / 2 a( d+1 > Ad v-, s j+l iji kb)[ii> / /ITTTA g sinh (^) - j cos g cosh(£%) 

J (o,o) - ^1/2 y J^- 1 ) e 7 ^^-( rf+ D/2^v / y + 1) y2 + cos2(/) 

[32l eqn.(3.28)], one concludes 



J = 



f(W) %M-%(«) 



,COS( 



— I 



a7 J(0 ' 0) 



Substituting into (|B.23|) one gets the recurrence relation: 



J (i -i) - J2 J (0 -1) + 



e 



icoscj) 



j, 



(0,0) 



The property is hold for all indices (r ^ —1), thus we can generalize: 

J(i,t) = ^J{o,t) + -£§ e^ cos ' l 'K d /2 +T+ i(z') - Kd/2+r+i{z) 
Splitting it onto z— and z'— parts gives 



20S(j> i <9J(oi-) 



(cos< 



M P M A dj> ( s ic°s<f> T [z] 

J (1,t) - £2 J (0,t) ^2"- fV d/2+r+U2J - J ( ,r+1) 



jW] _P_j[A .^djt****, (JS ic0i 

J (l,r) - £2 J (0,r) + £2 R A d/2+T+l(Z ) £ 



^ _ icosct) [z '] 

J (0,r+1) 



Differentiating (|B.3f I) over cos (f>, one gets 



J(2,r) - ^J(1,t) + r A (i /2 +r+ i(2: J - ^(O.r + l) £ J (l,r + 1) 



Substituting (|B~3T) 



^(2,t) =|4 >/(o,r) + ~j^-e lxiicos ' t 'K d/2+T+1 (( n ) 



1 2i/3 cos</> 



~. -Ao.r-! I) + 1je i4COS0 ^ rf / 2+T + l(2:')' 



(B.29) 



(B.30) 



(B.31) 



(B.32) 



(B.33) 



f ik d COS(j) „ ix £ cos4>T> , . 

|2 J (0,T+1) £3 e ^rf/2+T+2(Cnj 



£ 5 



— ^(0,r+2) • 



(B.34) 
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Splitting, 

j[z] _P 2 T [z] A d f3 ~ 2i(3 cos</> H 1 [z] iAdCosc/)-- cos 2 <j) M 

J (2 : r) -|4 J (0,r) ~ A d/2+r + l (Z J |g J (0,r+1) £2 J (0,r + 1) H £3 A d/2+r+2^J £2 _J (0,r+2) 



,r) ^^4 J (0,r) + A d— e - ft -d/2+T+lV* J £3 J (0,r + 1) £S 

^ e - t M/2+r+2l 2 J | ' 



— e^K d/2+T+2 (z>) - -^J^ T+2) ■ (B.35) 



Substituting (jR32l) into pOUj) gives the full expansion of vectorial integral in terms of Macdonald functions. 

Appendix C: Integration over frequencies and angles 

(a) Here are computed the angular integrals of the general form 

| g-/ (i-?«V ' M <C1) 

o 

with integers m, n, which were encountered in the text. For 2m > n + 1 one finds in the leading order [3ll | 

T/n f 1 V 2 7 x V" 2 / 2|7i-ll-l /pi <^ 

r(m) 



For 2m < n + 1 one obtains 



K = r \ 2 ; I, 2 , 1 ■ (C3) 

1 (n — m + lj 



In the case 2m = n + 1 the behavior of the integral is logarithmic. The only one of those needed here is the: 

4(ln2 7 -l) + 0( 7 - 2 ). (C.4) 



, 2 
V 3 — 

2 — -i 
V 6 



i l + v o 
In 2?7 

1 — v 



(b) Integration over the remaining angles is performed next. Since, with our choice of the coordinate system 
only one extra angle tp enters the amplitudes, only the following formula is needed: 



sD _ 3 1 ** v \»dn = r{D - l+N) r [— ) ■ (C5) 

(c) Finally, computation of the integrals over the frequency or over the impact parameter involving two Mac- 
donald functions of the same argument is performed using the formula (44[: 

K^cz)K v {cz)z a - x dz = K 2 > K 2 : K 2 jM 2 ' ■ (C.6) 



(i 



c a r(a) 



Appendix D: Asymptotic behavior of retarded fields 

In this Appendix we derive the asymptotic behaviors of the retarded solutions of the gravitational wave 
equation in arbitrary dimensions and of their space-time derivatives. 

The derivation is based on the well known formula [45^ for the retarded Green's function of the D— dimensional 
d'Alembert operator 

[^^)^ /2 - 2) (^), ^=even; 
with S^ix 2 ) = 6{x 2 ) and x 2 = {x ) 2 ~ {x 1 ) 2 - ... - (x ' 1 ) 2 . 
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First, suppose we have one point particle, moving along the trajectory T, parametrized by z (t) and 
properly normilized (i 2 = 1). The corresponding source of gravitational the field on flat background is 
T MN {x) = mj z M (t)z n (t) S(x - z(t))cLt. The solution of the d'Alembert equation nip*™ = -x D T MN is 
the convolution ip MN (x) = G rc t * T MN . Let us denote by f the retardation point on the smooth particle- 
worldline V corresponding to a given observation point x M = (t,x). Given x M and the time-like trajectory 
z m (t), t is the unique solution (with x° > z°(f)) of the equation (x — z(t)) 2 = 0, and specifies the intersection 
of r with the past light-cone of the observation point x M . Using hats to denote all corresponding kinematical 
quantities (e.g. z = z(j )), the Lienard-Wiechert solution becomes 27 

^MN,.s "Dm ( 1 d \^ l M l N 

where 



SF^T^Sf.' — (D - 2) 



1 d 

P&) = - 2 Qf( x ~ 2 ) 2 = k ■ ( x ~ 5 ) > • 

Note that in the co-moving at the retardation moment f Lorentz frame, the quantity p{x) coincides with the 
spatial distance |x — z| between the observation and the retardation points. 

Now, introduce the null vector c = (x M — z M )/p. It represents the properly normalized vector from the 
retardation point towards the observation point. Introduce in addition the vector n = c — i, which satisfies 
n 2 = — 1 and n ■ z — 0, as a consequence of z ■ c = 1 and z 2 = 1. Thus, c M is the sum of the orthogonal unit 
vectors: z , which specifies a time direction, and n M which is space- like and determines a spatial direction, 
i.e. a point on the sphere S D ~ 2 , which is the "basis" of a null-cone with axis along z M and apex at z M . Thus, 
for an arbitrary point x M on this light cone, p is the distance of the apex to the center of the sphere which 
contains x M , and is also equal to the radius of that sphere. Thus, going to the wave-zone of radiation emitted 
by the accelerated particle corresponds to considering the limit p —¥ +oo keeping fixed the quantities n, z, z, 
etc, which are related only to direction and refer to the apex z of the light-cone. 

Now with the help of the relations [46| 

9mt = c m P= {zc)p - 1 

9m P = z M + pen d M c N = - ($ M - z M c N - c M z N - pc M c N ^j (D.3) 

9m zn — cm 'zn , and similar ones for higher derivatives, 

and consecutive-differentiation rule (|D.2|) . one can check explicitly that in even dimensions the asymptotic 
behavior of ip MN and its derivatives is 

jMN lMN 

• I'd" = + - + ^Hf. with ^nld , ... . <ewt - tensors depending on c, z, 2, ... 

jUMNP \nMNP 

• i(>d N ' P = ^5jtr + - + with *Rld P > - ,*Newf " tensors depending on c, 1,1, ... . The same 

asymptotic behavior holds for higher order space-time derivatives, but not for Uip MN , for which one 
obtains instead 



• Dip 



jMN lMN 
MN _ yRad , , ^Ncwt 
p D/2 pD-X 



In odd dimensions, on the other hand, the solutions of d'Alembert equations are not known in closed form 
and one has to rely on approximate asymptotic expansions. In particular, one can use the expansion (proven 
for scalar and vector fields in [13] and conjectured for the tensor field ip MN ) 

j.MN 

^ D IN -% L T+°^ D/2 )^ (°-4) 



The square root derivative, which appears in odd-dimensions, is defined as the convolution 
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where ip^J^ again depends on c and on the retarded quantities z, z etc, but is expressed as an integral along 
the particle trajectory from r' — — oo to t' — f. Taking into account that c enters into ipnJd either in the 
product (z — z(t')) ■ c or with a free index, one can apply (|D.3I) to the integral and, using c 2 = to obtain the 
same behavior as in even dimensions, namely 28 : 



lMN 



Vp^=0(p- D ^) and U^ N =0{p- D ' 2 ). (D.5) 

For full gravity these properties do not hold exactly due to non-linearities, but to the order discussed here we 
have effective linear theory with linear operator □ and sources 1 T, 1 T I and S( l h, 1 h'), respectively. The first 
two are local currents and satisfy the assumptions made above. S( 1 h), on the other hand, is not a current 
generated by the point-like source, but, expressed explicitly with the help of (|A.5|) . and taking into account 
the field equation f|2 . 23[) . one can easily check that this part of □ 2 ipMN falls-off much faster at large distances, 
namely as the product of two Coulomb fields, 

2 ^mn(S) = -k d Smn ~ 1 h 1 ti = 0{p- {2D - 6) ). (D.6) 

Thus, when substituted to replace □ 2 4>mn into the integrals 

J 2 h PQ U 2 Tp M pda Q and J ' 2 hU 2 ^ MP da p 
c c 

encountered in Section 2, the resulting integrals vanish as a result of the asymptotic behaviors (|D.5|) and (|D.6|) . 
In other words, Smn is significant only in the space-time volume in which both h and h' are important, i.e. 
close to the moment of collision. This justifies the use of the asymptotic expansions presented above. 
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